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ABSTRACT 


Rainfall  and  othar  variables  with  similarly  skewed  distributions 
ara  hard  to  characterise  climatically  due  to  their  extreme  variabil¬ 
ity.  The  arithmetic  seen,  enarally  used,  ia  greatly  influenced  by 
extreme  values.  For  rainfall  data  Aron  219  station*  located  in  the 
western  United  States,  the  Median  was  found  to  be  a  nore  representa¬ 
tive  value,  and  somewhat  better  than  the  naan  for  predicting  future 
rainfall  aaounts.  Scsm  Monthly  precipitation  fr<  tumicy  distributions 
are  so  greatly  skewed  that  values  snaller  than  the  naan  occur  90Jt  of 
th*  tine.  Because  any  elngle  Measure  of  central  tendency  is  incon¬ 
clusive,  Measures  of  absolute  and  relative  variability  are  sumarised. 
Naps  of  percentage  occurrence  of  the  naan,  ratio  between  Median  and 
nean,  coefficient  of  variation  (CV),  and  relative  variability  (Vq)  are 
presented  for  the  nid-season  Months— Jan,  Apr,  Jul,  and  Oct,  ' 
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CHAPTER  I 


INTRODUCTION 

The  part  played  by  numerical  material  in  geography  has 
been  progressively  widened  in  recent  years,  but  the  methods 
of  generalizing  therefrom  have  remained  substantially  un¬ 
changed.  The  arithmetic  average,  or  mean,  is  still  the 
most  widely  used  value  to  characterize  a  set  of  observa¬ 
tions.  It  is  the  purpose  of  this  thesis  to  examine  the 
suitability  of  the  mean  as  a  measure  of  central  tendency  in 
non-symmetrical  distributions,  and  the  possible  advantages 
of  the  median  over  the  mean. 

The  introductory  chapter  is  concerned  with  the  reasons 
for  the  acceptance  of  the  mean,  and  with  previous  studies 
dealing  with  the  median  and  mean.  Chapter  II,  NClimatic 
Normals,"  presents  a  detailed  study  of  one  instance  where 
the  mean  is  being  used  to  represent  non-normally  distributed 
data.  Chapter  III,  "Measures  of  Dispersion,"  points  out 
the  limitations  of  any  single  value  in  portraying  the  dis¬ 
tribution  of  a  set  of  observations.  The  final  chapter. 
Chapter  IV,  "Climatic  Prediction,"  demonstrates  the  appli¬ 
cation  of  knowledge  gained  fro.n  previous  chapters  in  the 
solution  of  a  practical  problem. 

1.1  Limitation  of  Study 

Although  there  has  been  an  increasing  use  of  numerical 
analysis  in  climatology,  the  statistics  being  used  to  ana- 
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lyze  the  masses  of  data  have  gone  unchanged.  The  statistics 
presently  in  use  were  developed  for  use  with  normal  or  near 
normal  samples.  Chief  among  these  are  the  two  most  commonly 
used  statistics,  the  mean  and  the  standard  deviation. 

Unfortunately,  most  climatic  variables  have  distribu¬ 
tions  other  than  normal.  Singly  or  doubly  bounded  variables 
have  skewed  distributions  markedly  different  from  the  sym¬ 
metry  of  the  normal  distribution  (Section  2.5).  Precipita¬ 
tion,  wind  speed,  insolation,  and  visibility  are  a  few  of 
the  climatic  variables  for  which  the  majority  of  observa¬ 
tions  occur  at  the  upper  and/or  lower  limits  of  their  dis¬ 
tributions.  In  such  skewed  distributions  the  mean  and  stan¬ 
dard  deviation  do  not  have  all  the  useful  properties  as  in 
the  normal  distribution.  Perhaps  the  median,  and  measures 
of  dispersion  about  it,  may  be  more  useful. 

Although  the  discussion  will  be  limited  to  the  use  of 
the  mean  and  median  in  summarizing  precipitation  data,  the 
results  should  be  equally  valid  for  other  non-negative  vari¬ 
ables  having  substantial  frequencies  near  the  lower  bound, 
and  thus  having  similarly  skewed  frequency  distributions. 
Wind  speed,  income,  and  agricultural,  industrial,  or  mineral 
production  by  counties,  are  good  examples  of  such  variables. 

1.2  Use  of  Mean  and  Median 

The  arithmetic  mean  is  often  declared  to  be  the  bast 
estimate  of  central  tendency,  and  hence  is  the  value  most 
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generally  used  to  characterize  entire  sets  of  observations. 
In  a  symmetrical  distribution  such  as  the  normal,  the  mean 
is  also  the  median  or  50#  probability  value.  Thus  many 
people  implicitly  regard  any  mean  as  having  a  50#  probabil¬ 
ity  level  of  exceedance.  Unfortunately,  in  skewed  distri¬ 
butions,  which  are  more  common  for  geographic  data,  the 
mean  is  not  the  median,  and  only  the  median  has  a  50#  pro¬ 
bability  level  of  occurrence. 

The  greater  the  skewness  of  the  distribution,  the 
poorer  the  mean  becomes  as  an  estimate  of  the  middle  value, 
the  median.  While  the  mean  may  be  easier  to  work  with  in 
many  instances,  the  median  or  50#  probability  value  is  most 
often  desired  and  should  be  used  whenever  possible  to  avoid 
any  misconceptions  regarding  frequency  of  occurrence.  In 
many  applications  the  mean  has  been  uced,  and  is  continuing 
to  be  used,  largely  as  an  approximation  of  the  median. 

1.3  Early  Studies 

The  unsuitability  of  the  mean  a  a  measure  of  central 
tendency  in  climatological  research  has  long  been  realized, 
but  it  was  not  until  1933  that  a  serious  attempt  was  made 
to  replace  the  mean  with  the  median  to  summarize  a  "varied 
meteorological  record."  It  was  in  this  year  that  the  Bri¬ 
tish  geographer,  Percy  Robert  Crowe,  introduced  the  use  of 
the  median  and  quurtile  deviation  (half  the  difference  bet¬ 
ween  the  upper  and  lower  quartiles)  in  a  ^tudy  cf  European 
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rainfall  entitled  "Analysis  of  Rainfall  Probability."  This 
study  created  much  interest  and  its  influence  can  be  seen  in 
several  succeeding  studies.  The  most  notable  of  these  is  a 
study  of  Indian  rainfalls  by  Matthews  (1936).  The  median 
and  quartile  deviation  were  again  used  by  Crowe  (1936)  to 
study  the  rainfall  regime  of  the  Western  Plains  of  the 
United  States,  and  also  by  Lackey  (193V)  in  constructing 
annual-variability  .naps  of  the  Great  Plains.  Earlier  stud¬ 
ies  by  Lackey  (1935,  1936)  also  show  the  influence  of  Crowe. 

Other  researchers  not  influenced  by  Crowe  were  also 
coming  to  the  conclusion  that  the  mean  should  be  replaced 
by  the  median  when  dealing  with  a  highly  variable  distri¬ 
bution.  Gisborne  (1935)  used  monthly  precipitation  records 
of  Spokane,  Washington  to  show  that  the  mean  is  a  poor 
value  to  use  as  a  "climatic  normal."  He  found  that  the 
mean  was  never  reached  or  exceeded  more  than  46%  of  the 
time  in  any  month  and  that  in  June  the  mean  occurred  only 
26 %  of  the  time.  If  the  mean  were  used  as  the  "normal," 

74  years  out  of  100  had  below  normal  precipitation  in  June. 

A  similar  study  was  undertaken  by  Mindiing  (1940)  using  60 
years  of  monthly  precipitation  records  for  each  of  14  sta¬ 
tions  distributed  throughout  the  United  States.  His  con¬ 
clusions  were  the  same  as  Gisborne's:  the  mean  is  an  unsat¬ 
isfactory  value  to  use  as  a  "normal,"  and  that  the  median 
would  be  a  better  value  to  use. 

This  interest  in  the  proper  statistic  for  use  in 
representing  a  non-symmetrical  distribution,  first  seen  in 
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the  studies  by  Crowe,  reached  a  peak  June  20-21,  1940.  At 
this  time  a  resolution  recommending  that  "the  expression  of 
normals  of  precipitation  in  future  hydrologic  studies  be 
defined  by  the  median  instead  of  the  arithmetical  average" 
was  presented  to  the  Section  of  Hydrology  of  the  American 
Geophysical  Union,  meeting  jointly  with  the  Western  Inter¬ 
state  Snow-Survey  Conference  in  Seattle.  A  committee  headed 
by  P.E.  Church  (1941)  had  weighed  the  advantages  and  dis¬ 
advantages  of  each  (Chart  I)  and  concluded  that  "in  the 
future,  at  least  for  hydrologic  studies,  the  expression  of 
normals  of  precipitation  be  defined  by  the  median." 

Unfortunately  the  recommendations  of  this  committee 
were  lost  with  the  on-coming  of  World  War  II,  as  was  the 
interest  stimulated  by  Crowe  and  others,  and  the  mean  con¬ 
tinues  to  be  used  in  hydrologic  and  climatological  studies. 

1.4  Recent  Studies 

Very  little  discussion  of  this  topic  has  appeared 
since  1941.  Land^berg  (1947),  declaring  that  the  mean  does 
not  represent  the  usual  condition  even  fcr  monthly  tempera¬ 
tures,  condemned  the  use  of  the  term  "normal"  for  the  arith¬ 
metic  mean.  In  a  study  of  rainfall  on  Oahu,  Hawaii,  Lands- 
berg  continued,  saying  that: 

The  median  is  the  statistic  that  preserves 
one  6f  the  most  important  properties  of  a  "nor¬ 
mal."  It  represents  the  center  of  the  distribu¬ 
tion  and  half  of  the  observations  are  higher, 
the  other  half  lower  -than  this  value.  (19bl,  p.  9) 
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In  the  study  of  the  measurement  of  water  content  of  a 
snow  course,  Court  (1958)  found  that  a  water  equivalent 
figure  sufficiently  accurate  for  the  practical  use  to  be 
made  of  it  can  be  obtained  from  the  median  of  a  small  num¬ 
ber  of  snow  course  points.  He  indicates  that  the  combina¬ 
tion  of  a  larger  number  of  snow  water  content  measurements 
into  a  single  mean  is  "a  gross  waste  of  information  from 
the  statistician's  standpoint." 

More  recently  Joos  (1964),  after  a  study  of  the  vari¬ 
ability  of  weekly  rainfall,  suggested  that  "a  more  meaning¬ 
ful  ’normal*  for  weekly  precipitation  might  be  based  on  the 
50#  value  rather  than  on  a  long  term  mean.”  Kangieser 
(1966)  finds  that  in  some  parts  of  Arizona  only  one  year  in 
ten  has  monthly  precipitation  as  great  as  the  mean  value. 
And  Bennett  (1967)  shows  the  median  to  be  a  better  statis¬ 
tic  than  the  mean  for  the  study  of  insolation  which,  un¬ 
like  precipitation,  has  a  distribution  in  which  many  obcer- 
vations  occur  near  the  upper  bound. 

1.5  The  Mean  as  a  Poor  Estimate  of  Central  T  nd*-ncy 

The  mean  is  too  strongly  influenced  by  the  extremes  in 
a  series  of  values.  Ten  months  with  values  one  unit  below 
the  mean  are  required  to  neutralize  the  effect  of  one  month 
with  a  value  ten  units  above  the  mean.  In  the  case  of  cli¬ 
mate,  it  is  the  time  period  during  which  people  have  to 
live  and  work,  the  month  of  the  year,  that  is  significant. 
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The  frequency  of  departures  from  the  normal  is  important, 
but  the  actual  extent  of  departure  of  each  instance  is 
secondary. 

Yet,  it  is  the  mean  which  is  stated  in  answer  to  a 
question  concerning  the  amount  of  rainfall  a  place  receives. 
At  the  present  time,  maps  of  climatic  normals  are  based  on 
the  arithmetic  mean  such  as  those  which  appear  in  The  Na¬ 
tional  Atlas  of  the  United  States,  published  in  1964.  The 
supposed  purpose  of  these  "normal  total  precipitation" 
maps  is  tc  give  the  user  an  idea  of  the  most  likely  amount 
of  precipitation  to  be  received  at  a  particular  place. 

Since  the  arithmetic  mean  can  satisfy  the  concept  of  normal 
only  when  the  data  are  symmetrically  distributed,  these 
"normal  total  precipitation"  maps  are  deceptive  especially 
in  drier  areas  of  this  country. 

It  would  seem  that  the  median  would  be  a  much  more 
logical  measure  of  central  tendency.  Since  it  is  the  mid¬ 
dle  value,  it  is  easily  ascertained  when  the  figures  are 
arranged  in  ascending  or  descending  order  of  magnitude  and 
can  thus  be  located  by  a  process  no  more  complicated  than 
simple  inspection  or,  at  most,  by  plotting  the  values  con¬ 
cerned  on  some  form  of  linear  graph.  With  the  use  of  the 
high  speed  computers,  even  large  amounts  of  data  pose  no 
problem  in  the  finding  of  the  median.  Its  reality  consists 
of  the  fact  that,  since  half  of  the  records  exceed  it  and 
half  of  them  are  less,  it  represents  the  usual  or  typical 
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value.  It  is  quite  independent  of  a  few  very  exceptional 
values,  since  each  value  has  an  equal  pull,  whatever  its 
size.  Because  of  this,  the  median  is  less  sensitive  to 
errors  which  might  be  present  in  the  data  being  examined. 

1.6  Normals 

At  the  present,  one  of  the  primary  methods  used  to 
portray  climatic  data  is  through  the  use  of  a  'climatic 
normal.'  Climatic  normals  supposedly  describe  the  climate 
of  a  place  or  region  and  are  often  used  to  estimate  future 
climatic  conditions.  They  are  calculated  by  averaging 
climatic  events  over  a  number  of  years,  usually  30.  Al¬ 
though  some  objections  have  been  voiced  against  this  method 
of  determining  the  normal,  it  is  accepted  as  being  the  best 
possible  method  by  the  majority  of  the  people. 

As  pointed  out  by  Gisborne  in  1935,  the  concept  of 
normal  is  well  established  in  both  the  technical  and  the 
lay  mind  as  that  which  is  common,  natural,  ordinary,  regu¬ 
lar,  typical,  and  usual.  It  must  be  determined  just  how 
well  the  mean  satisfies  this  definition. 
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CHAPTER  II 


CLIMATIC  NORMALS 

2.1  Data  and  Study  Area 

Used  in  this  study  are  monthly  precipitation  values  at 
219  stations  in  the  United  States,  all  west  of  the  Missis¬ 
sippi  River  (Map  1).  Data  were  taken  from  publications  of 
the  United  States  Weather  Bureau,  supplemented  by  data  pub¬ 
lished  by  various  states.  Stations  were  selected  to  insure 
homogeneous  data:  no  significant  change  in  station  posi¬ 
tion,  elevation,  or  environment  during  the  observation 
period.  Although  some  stations  had  precipitation  records 
for  110  years  of  continuous  observations,  only  data  from 
the  period  1931  to  1960  were  used.  Thus  all  analyses  and 
maps  in  this  report  are  comparable,  and  are  for  the  period 
specified  by  the  World  Meteorological  Organization  and  used 
by  the  United  States  Weather  Bureau  to  determine  climatic 
normals . 

Table  1  (Appendix  1)  is  a  list  of  the  monthly  preci¬ 
pitation  stations  used  in  this  study.  The  219  stations 
represent  an  average  density  of  one  station  for  every  9,789 
square  miles.  While  this  falls  short  of  the  one  station 
per  1,029  square  miles  coverage  for  rainfall  by  the  Nation¬ 
al  Atlas  of  the  United  States ,  it  compares  favorably  with 
the  density  of  stations  of  the  Local  Climatological  Data 
published  by  the  United  States  Department  of  Jomn.erce  Wea- 
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ther  Bureau  (approximately  one  station  per  10,500  square 
miles).  Stations  are  listed  by  their  U.S.  Weather  Bureau 
identification  number,  with  name,  coordinates,  and  number 
of  years  of  continuous  record.  For  convenience  in  compu¬ 
tation  and  mapping,  stations  are  also  numbered  serially 
1  to  219. 

Due  to  the  excessive  time  and  work  involved  in  mapping 
the  results  of  all  months,  only  the  mid-season  months  of 
January,  April,  July  and  October  are  included  in  this  study. 
The  use  of  these  months  should  bring  out  any  seasonal  vari¬ 
ations  that  may  exist,  and  should  be  a  large  enough  sample 
to  work  with.  It  is  important  to  work  with  more  than  one 
or  two  months  since  the  possibility  does  exist  that  they 
may  not  be  representative  of  the  majority  of  months.  Such 
an  error  was  made  by  McClean  (1956)  in  assuming  that  con¬ 
clusions  reached  for  one  month  would  apply  to  all  months. 
McClean' s  contention  that  the  period  1881-1915  is  a  poor 
one  to  use  as  the  standard  for  reference  to  British  rain¬ 
fall  was  disproved  by  Giasspoole  (1956). 

2.2  Procedure 

If  the  Limatic  normal  were  a  representative  value,  it 
would  be  expected  to  be  reached  in  half  of  the  years  of 
record.  The  percentages  of  the  years  (1931-1960)  in  which 
the  "normal"  value  was  reached  or  exceeded  were  calculated 
and  are  shown  on  Map  2-5.  These  maps  show  that  the  mean 
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fails  to  represent  the  central  value  of  the  distribution 
from  which  it  is  derived,  especially  for  the  Southwestern 
portion  of  the  country.  The  inadequacy  of  the  mean  as  an 
expression  of  the  "normal"  for  monthly  precipitation  is 
shown. 

.?.3  Skewness 

As  stated  in  the  introduction,  in  a  normal  distribu¬ 
tion,  or  any  other  symmetrical  distribution,  the  mean  is 
equal  to  the  median  and  under  such  circumstances  the  mean 
is  a  perfectly  valid  value  to  use.  But  when  the  distribu¬ 
tion  is  skewed  (asymmetrical),  the  mean  is  much  less  repre¬ 
sentative.  A  distribution  is  symmetrical  when  the  median 
equals  the  mean;  positively  skewed  when  the  mean  is  larger 
than  the  median;  and,  negatively  skewed  when  the  mean  is 
smaller  than  the  median. 

There  are  several  ways  for  measuring  skewness.  One 
of  these  is  the  "Pearsonian  coefficient  of  skewness"  given 
by  the  formula 

3(mean  -  median)/standard  deviation 
Another  measure,  <s  (alpha-three),  is  defined  in  terms  of 
the  second  and  third  moments,  m2  and  m3,  about  the  mean 
as 

«» •  %/(■«>* 

In  this  study  a  simple  measure  of  skewness  is  given  by  the 
percentage  occurrence  of  the  mean  as  portrayed  in  Maps  2-5. 
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The  further  the  mean  departs  from  the  50?$  probability  value, 
the  greater  the  skewness  and  the  greater  the  advantage  in 
using  the  median  to  best  represent  what  is  expected  to 
happen  again. 

2.4  Analysis  of  Maps 

Before  discussing  the  maps  any  further,  it  should  be 
pointed  out  that  isolines  are  only  as  good  as  the  data  from 
which  they  are  derived.  The  patterns  shown  on  the  maps  in 
this  study  are  necessarily  general  due  to  the  number  of 
stations  used.  If  further  data  points  were  available,  pat¬ 
terns  would  undoubtedly  be  changed.  But  for  the  purpose  of 
the  present  study  the  number  of  stations  is  adequate.  If 
further  detail  is  desired,  studies  of  individual  states 
should  be  made,  such  as  that  of  Arizona  by  Kangieser  (1966). 
It  would  have  been  of  interest  to  compare  his  maps  with  the 
ones  in  this  study,  but  unfortunately  they  are  of  different 
months . 

All  maps  show  low  values  to  occur  with  regularity  in 
the  Southwest,  while  high  values  occur  in  the  Northwest. 

The  mean  is  reached  in  as  many  as  50?$  of  the  years  only 
for  January  in  the  northwestern  section  of  Washington.  If 
the  mean  were  a  good  value  to  use  as  a  climatic  "normal," 
it  should  have  been  reached  in  half  the  years  in  all  parts 
of  the  countiy.  But  it  is  reached  in  less  than  a  quarter 
of  the  years  in  many  areas  of  the  country,  and  in  fewer 
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than  one  year  in  ten  for  the  Central  Valley  of  California 
in  July. 

Seasonal  variation  in  the  occurrence  of  the  mean  can 
best  be  seen  in  central  California.  Here  the  values  range 
from  less  than  10#  in  July  to  over  4b#  in  January. 

For  all  months  values  from  35#  to  45#  are  most  common. 
There  can  be  seen  no  general  longitudinal  or  latitudinal 
gradation  of  values,  and  the  only  consistencies  are  the  low 
values  of  the  southern  California  desert  and  the  high 
values  of  northwestern  Washington. 

2.5  Conversion  of  Existing  Maps 

If  a  distribution  proves  to  be  skewed,  the  use  of  the 
mean  as  the  most  representative  value  of  the  distribution 
will  be  misleading.  Under  such  circumstances  the  median 
should  be  used.  However,  considerable  work  went  into  the 
computation  of  existing  means,  and  it  is  doubtful  that  any¬ 
one  would  be  willing  to  discard  these  even  if  their  valid¬ 
ity  is  in  doubt.  This  problem  would  be  solved  if  a  stable 
relationship  is  found  to  exist  between  median  and  mean. 

With  this  aim  in  mind,  the  ratio  between  median  monthly 
and  mean  monthly  values  was  computed  for  each  station 
(Maps  6-10).  If  stable,  these  ratios  can  be  used  as  cor¬ 
rection  factors  to  convert  existing  means  to  medians  (see 
Appendix)  . 
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RATIO  OF  MEDIAN  TO  MEAN 
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The  ratios  between  the  median  and  mean  can  also  be 
used  as  excellent  measures  of  skewness.  A  ratio  smaller 
than  1  indicates  a  positively  skewed  distribution,  while  a 
ratio  larger  than  1  indicates  a  negatively  skewed  distribu¬ 
tion.  The  greater  the  departure  from  1,  the  greater  the 
degree  of  skewness. 

As  would  be  expected  these  maps  compare  favorably  with 
maps  2-5  with  low  values  occurring  in  the  Southwest  and 
high  values  in  the  Northwest. 

Several  stations  in  California  had  no  rain  at  all 
during  the  30  years  studied  for  the  month  of  July  and  had 
to  be  omitted  or  overlooked  in  the  drawing  of  isolines, 
since  in  this  case  the  ratios  would  have  been  infinite. 
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CHAPTER  III 


MEASURES  OF  DISPERSION 

3.1  Absolute  vs.  Relative  Measures  of  Variability 

Both  the  mean  and  the  median  have  their  limitations 
in  that  they  are  measures  of  central  tendency  only.  No 
single  value,  however  determined,  can  hope  to  give  a  com¬ 
plete  descriptive  summary.  Crowe  (1933)  pointed  out  that 
we  require  not  only  an  index  of  the  "normal  sequence  of 
events,"  but  also  of  the  "probable  frequency  and  extent  of 
variations  from  the  normal." 


3.2  Absolute  Measures  of  Variability 

The  standard  deviation  is  used  to  estimate  the  degree 
of  variability  about  the  mean  and  defined  by  the  equation 


S.d. 
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where  n  is  the  number  of  values,  is  the  individual  value, 
and  x  the  mean  of  the  n  values.  But  the  standard  deviation 
is  not  an  appropriate  t^rm  to  be  u^ed  with  the  median: 
while  the  sum  of  the  squared  deviations  is  least  when  com¬ 
puted  about  the  mean,  the  sum  of  the  absolute  deviations  is 


at  a  minimum  when  computed  about  the  median.  A  measure  of 
dispersion  used  to  measure  variability  about  the  mean  or 
median  ic  the  mean  deviation,  also  known  ac  the  average 
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variability  (AV) .  The  mean  deviation  is  defined  by  the 
equation 

1 

m .  d .  =  n 

1  =  1 

where  z  may  be  either  the  mean  or  median.  The  semi-inter¬ 
quartile  range  or  quartile  deviation  q,  defined  by  the  equa¬ 
tion 

q  =  (Q3  -  Qx)/2 

where  and  Q3  are  the  first  and  third  quartiles,  is  used 
to  estimate  the  degree  of  variability  about  the  median. 

Other,  but  poor  measures  of  variability,  use  only  the 
highest  and  lowest  recorded  values,  M  and  m.  The  range, 

M  -  m,  is  often  used.  Heilman  (1909)  defined  the  "ratio  of 
variation"  as  the  quotient  M'm.  These  measures  are  not 
based  on  all  observations,  and  the  ratio  of  variation  can¬ 
not  be  used  for  desert  rainfall:  M/m  would  become  infinite 
for  many  stations. 

3.3  Relative  Measures  of  Variability 

Absolute  measures  of  variability  are  useful  in  helping 
to  understand  a  particular  set  of  observations,  but  do  not 
give  a  complete  picture  of  the  variability.  Thus  they  are 
of  little  value  when  comparisons  between  ooservations  at 
several  different  localities  are  desired.  This  is  due  to 
the  fact  that  variability  generally  increases  as  the  values 
of  the  observations  become  larger.  Therefore,  a  map  showing 
the  dispersion  of  observations  about  the  mean  or  median  at 
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various  stations  using  any  of  the  measures  of  variability 
listed  above,  would  give  the  same  picture  as  a  map  of  the 
means  or  medians  of  the  observations  themselves. 

For  this  reason,  relative  measures  of  variability  must 
be  considered  in  order  to  derive  comparable  figures.  The 
four  measures  of  relative  variability  corresponding  to  the 
first  four  measures  of  variability  mentioned  above  are 
given  in  Table  2. 

Of  these  measures,  only  CV  and  Vr  have  been  studied  to 
any  great  extent.  In  studying  the  relative  variability  (Vr) 
of  precipitation,  Conrad  (1941)  found  that  a  very  strong 
mathematical  dependence  of  the  value  of  Vr  on  the  yearly 
sum  occurs  when  these  sums  are  below  1000mm  (39.4  inches). 
Because  of  this  he  concluded  that  the  relative  variability 
(\/r)  cannot  be  used  to  compare  the  variability  of  precipi¬ 
tation  in  a  locality  with  a  small  annual  fall  with  that  in 
another  locality  where  the  precipitation  is  large. 

In  a  study  of  all  four  measures,  Longley  (1952)  lound 
the  coefficient  of  variation  the  most  satisfactory  measure 
for  comparing  variability  of  precipitation  between  dif¬ 
ferent  stations.  It  was  found  to  have  small  errors  through 
sampling,  in  that  the  values  of  CV  are  comparable  even  when 
the  mean  precipitation  values  are  quite  different. 

Longley  also  found  that  variability  tends  to  be  greater 
where  the  precipitation  is  least,  but  the  relationship  was 
not  close.  In  a  study  of  precipitation  data  at  34  stations 
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TABLE  2.— RELATIVE  MEASURES  OF  VARIABILITY 


The  Coefficient  of  Variation: 


The  Relative  Variability: 


and 


VQ 


100  (Q3  -  QJ 


and,  also 


Qr  .  M_-_m 

z 


n  is  the  number  of  values 
is  the  individual  value 
R  is  the  mean 
T  is  the  median 
z  is  the  mean  or  median 
M  is  the  maximum  value 
m  is  the  minimum  value 


in  British  Columbia  and  Washington,  he  found  the  coefficient 
of  correlation  between  mean  precipitation  and  the  coeffi¬ 
cient  of  variability  to  be  -0.68  for  July,  -0.71  for  Decem¬ 
ber,  and  -0.48  for  the  annual  data. 

3.4  Variability  as  Related  to  Rainfall  Amounts 

Since  the  majority  of  stations  studied  in  this  report 
have  annual  rainfall  amounts  below  30  inches,  it  is  impor¬ 
tant  to  know  the  relation  between  Vr  (and  CV)  or  Vq  and  the 
mean  or  median  rainfall.  Conrad  states  that: 

".  .  .the  assumption  that  Vr  represents  a 
numerical  characteristic  of  variability,  unre¬ 
lated  to  the  arithmetic  mean,  has  been  proven 
fallacious  by  the  observations.  Therefore,  con¬ 
clusions  drawn  comparing  values  of  Vr  for  dif¬ 
ferent  places  in  the  vast  regions  where  the 
annual  precipitation  is  less  than  28  inches.  .  . 
are  inaccurate  and  misleading."  (in  Conrad  and 
Poliak,  1950,  p.  56) 

Longley's  study  indicates  that  this  may  not  be  the  case. 

An  attempt  will  here  be  made  to  clarify  this  question  and 
to  also  look  at  monthly  relationships  where  the  rainfall 
amounts  are  still  smaller. 

Figures  1,  2,  and  3  were  prepared  to  show  annual  end 
monthly  relations  between  CV  and  mean  rainfall,  and  Vq  and 
the  median  rainfall.  If  a  least  fit  line  were  fitted  to 
each  set  of  data,  a  curve  would  be  obtained  similar  to  the 
one  found  by  Conrad  when  examining  annual  rainfall.  The 
turning  point,  best  seen  in  July,  occurs  at  about  two  inches 
for  monthly  data  and  twenty  inches  for  annual  data. 
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Even  though  the  scattering  is  great,  at  first  glance 
it  would  seem  as  though  Conrad's  statement  was  true:  at 
lower  values  there  is  an  inverse  relation  between  the  mea¬ 
sure  of  variability  and  the  measure  of  central  tendency. 

But  these  figures,  as  those  by  other  researches,  overlook 
several  things.  When  performing  statistical  analysis  on 
data  we  cannot  take  samples  from  all  areas  (regions)  and 
treat  them  as  coming  from  one  'population'.  In  the  past, 
many  studies  have  been  made  to  derive  generalizations  about 
certain  climatological  variables.  In  these  studies  formu¬ 
las  were  derived  by  looking  at  data  from  many  stations 
scattered  throughout  the  world.  If  these  formulas,  or 
generalizations,  were  formulated  for  the  purpose  of  repre¬ 
senting  a  base,  and  the  regions  determined  by  how  they 
deviate  from  this  base,  this  is  a  valid  approach.  But,  if 
the  purpose  was  to  derive  a  generalization  from  which  pre¬ 
dictions  could  be  made,  the  study  was  invalid.  Generali¬ 
zations  formed  by  looking  at  data  from  dissimilar  areas 
can  be  applied  to  the  whole  mass  of  data,  but  rarely  to  any 
of  the  "subregions". 

Similar  applications  have  been  made  in  other  areas  of 
geography,  correlations,  regression  analysis,  and  other 
statistical  analyses  have  been  performed  on  samples  suppos¬ 
edly  from  the  same  'population'.  In  reality  these  samples 
have  come  from  several  different  'populations,'  as  is  the 
case  in  the  present  analysis. 


35 


In  the  majority  of  previous  studies  of  CV,  Vr,  and  Vq, 
data  from  widely  scattered  stations  in  regions  with  vastly 
different  rainfall  characteristics  were  treated  as  coming 
from  the  same  population  (or  region) .  The  results  of  these 
studies  supposedly  represented  the  behavior  of  these  mea¬ 
sures  of  relative  variability  for  all  areas.  That  this  is 
hardly  fact  can  be  seen  by  looking  at  the  behavior  of  these 
measures  separately  for  several-  different  areas.  Such  a 
study  is  presented  in  Figure  4  where  the  populations  or 
subregions  are  represented  by  A,  B,  C,  and  D.  The  bottom 
graph  in  Figure  4  shows  the  typical  curve  obtained  when 
data  from  many  different  populations  are  grouped  together. 
It  can  be  seen  that  a  generalization  arrived  at  by  consid¬ 
ering  populations  A,  B,  C,  and  D  grouped  together  would  not 
necessarily  hold  true  for  each  individual  population.  A 
large  inverse  relation  between  CV  and  the  mean  rainfall 
occurs  only  with  population  A  and  this  is  not  surprising. 
Any  measure  of  relative  variability  would  be  expected  to 
increase  rapidly  as  the  average  rainfall  approaches  zero, 
since  any  number  divided  by  zero  is  infinitely  large.  Al¬ 
though  only  CV  is  discussed  in  this  example,  the  same  rela¬ 
tionships  were  found  for  Vr  and  Vq  to  varying  degrees. 

3.5  Variability  Maps 

While  many  maps  of  mean  annual  and  mean  monthly  preci¬ 
pitation  for  the  United  States  exist,  few  maps  exist  cover- 
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ing  extended  regions  which  show  some  measure  of  dispersion 
of  the  annual  or  monthly  totals  around  their  means.  Two 
maps  were  uncovered  showing  the  annual  values  of  the  coef¬ 
ficient  of  variation;  the  first  by  Hazen  (1916),  and  the 
second  by  Hershfield  (1962)  (see  Map  11).  No  maps  were 
found  showing  monthly  rainfall  variability  although  the 
need  for  such  maps  was  realized  by  Hershfield.  Maps  12-15 
showing  monthly  values  of  the  coefficient  of  variation  were 
prepared  to  partially  meet  this  need  and  to  show  the  relia¬ 
bility  of  the  mean. 

No  maps  were  found  of  extended  regions  showing  the 
dispersion  of  values  about  the  median.  While  the  coeffi¬ 
cient  of  variation  is  a  suitable  term  to  use  in  measuring 
the  dispersion  around  the  mean,  a  good  index  of  variability 
around  the  median  is  obtained  by  •  -.pressing  the  quartile 
deviation  as  a  percentage  of  the  median.  Maps  16-19  were 
prepared  using  values  determined  in  thir  manner. 

3.o  A  Note  on  the  Coefficient  of  Variation 

The  coefficient  of  variation  (JV)  becomes  meaningless 
in  terms  of  normal  probability  theory  wh--n  it  become'  larger 
than  1.00  since  the  standard  deviation  oi  a  normally  distri- 
bited  variable  cannot  be  larger  than  the  mean.  But,  if  CV 
i’  looked  at  as  a  measure  of  dispersal  divided  by  a  measure 
of  central  tendency,  without  reference  to  th  type  ol  dis¬ 
tribution,  it  becomes  very  u.eful. 
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Hastings  (1965)  suggests  that  CV  can  be  used  for 
climatological  and  ecological  classifications,  and  as  an 
indicator  of  rainfall  probability.  If  seasons  are  ranked, 
the  one  with  the  lowest  coefficient  being  1  and  the  highest 
4;  the  rank-sequence  can  be  used  to  delineate  regions. 

Such  a  classification  reflects  accurately  tne  relative  va¬ 
lue  of  moisture  in  that  it  takes  into  account  not  only  sea¬ 
sonal  variability  and  amount  of  rainfall,  but  also  the 
reliability  of  th^  amount. 

The  coefficient  ot  variation  is  also  a  measure  of 
skewness.  Where  CV  exceeds  0.5  (50%)  skewing  sets  in. 

When  CV  exceeds  1.00  (100%)  non-normality  increases  rapidly. 
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CHAPTER  IV 


CLIMATIC  PREDICTION 


4.1  Normals 

In  previous  chapters  the  question  of  what  statistic 
best  represents  the  ’normal’  was  investigated.  In  this 
chapter  the  length  of  record  upon  which  to  base  the  normal 
will  be  examined. 

As  mentioned  earlier,  one  of  the  prime  uses  of  the 
normal  is  as  a  predictor  of  future  events .  At  present, 
climatic  normals  are  based  on  3C  yearr  of  record.  Many 
studies  (as  summarized  by  Court,  1967)  show  this  may  not  be 
the  optimum  length  of  record.  In  the  f  studies,  means  of 
varying  lengths  of  record  were  used  in  determining  the  op¬ 
timum  length  of  record.  No  attempt  har  y*.-t  been  found  to 
use  the  median  in  relation  to  this  question.  In  the  pre¬ 
sent  study,  the  median  o*  varying  lengths  of  record  will.  bo¬ 
used  as  the  predictor  of  future  eventr  ana  the  results  com¬ 
pared  with  those  obtained  through  the  use  of  the  mean. 


4.2  Procedures 

One  of  the  most  general 
mining  the  optimum  length  of 
averages  cf  the  squared  diff 
varying  numbers  of  observati 
dieted.  Using  the  notation 


y  u 

sed 

pro^edur 

,  i  f 

for  deter- 

rec 

ord 

i’ 

to  :a 

IlU 

late  the 

ren 

C'“S 

be' 

tween 

the 

means  of 

ns 

and 

tt| 

S  Valu 

being  pre- 

dop 

t«  d 

by 

Court 

(1 

Vo’)  thir 

49 


t 


extrapolation  variance,  3  ^m,  is  computed  as 


km 


1  V 

n-k-m+1  / 


n-k-m+1 


TTT 


k-1 


-i  * 


xi+krm-l  “  A  V” 

k  /  x ,  .  • 

F5  1+J 


(l) 


where  k  is  the  number  of  antecedent  observations  from  which 
a  mean  (or  median)  is  computed;  m  is  the  lag  between  the  k 
year  period  and  the  value  being  estimated;  and  n  the  number 
of  observations  x^  ordered  in  time  from  x^  to  xn. 

In  the  present  study  this  equation  was  modified  to 
give  D^,  the  average  of  the  absolute  differences  between 
the  medians  (U)  of  k  years  and  the  values  being  estimated, 
a  s 


n-k-m+1 


Jkm  = 


n-k-m+1 


V 


TTT 


xl+-k+m-l  "  xk,i 


(2) 


where  7^  j.  represents  the  median  of  k  years  beginning  with 
the  observation  x^. 

Since  Sj^  and  D^m  are  not  comparable  (the  first  being 
an  average  of  squared  differences  and  the  second  being  an 
average  of  absolute  differences),  a  third  value,  ,  the 
mean  prediction  error,  was  calculated  as 


n-k-m+1 


Q|tm  =  FIFSTT  y 

in 

As  Q^m  represents  tne  average 
ween  the  means  of  k  years  and 
it  is  directly  comparable  to 
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•  •  Appendix  II  contains  a  program  which  computes  the 
three  above  values  and  places  a  star  after  the  lowest  value. 
The  length  of  record  at  which  the  lowest  value  occurs  will 
be  hereafter  referred  to  as  k*.  The  program  also  computes 
the  percentage  difference  between  the  values  obtained  for 
S£m,  and  Q|<m  at  k*,  the  optimum  length  of  record,  and 

those  obtained  at  the  other  values  of  k. 

Although  the  program  computes  all  values  of  Sjcm  and 
Qj<m  for  k  years  from  1  to  50,  only  odd  values  of  k  from  3 
to  49  were  used  in  the  computation  of  Dj<m  in  order  to  keep 
the  program  relatively  simple.  For  this  reason,  and  in 
order  to  keep  everything  as  comparable  as  possible,  only 
the  values  of  s£m,  D^m*  and  Qkm  obtained  for  odd  values  of 
k  will  be  considered  in  most  of  the  following  discussions. 

4.3  Random  Numbers 

has  been  shown  by  Court  (1967)  to  decrease  as 
1  +  1/k,  and  the  graphs  appearing  in  his  report  bear  this 
line  as  a  guide  to  determine  the  significance  of  the  cal¬ 
culated  values  of  S^.  Since  is  no  longer  the  only  value 
being  examined,  it  would  be  of  interest  to  see  how  these 
lines  compare  with  those  for  ani  Dk  (m  no  longer  appears 
as  a  subscript  since  only  the  values  for  m~l  are  being  used.) 
To  do  this,  equations  (1),  (2),  and  (3)  were  used  to  com¬ 
pute  Si  Qk»  and  for  1000  random  normal  numbers  (mean 
zero,  variance  unity).  Values  were  also  obtained  by 
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TABLE  3  .--OPTIMUM  LEMQ3H  OF  RECORD  (k*)~ IRECIPITATION1 


STATION 

•JAN 

J 

riB 

_MAR 

.APR 

-MAX 

.JUN 

n 

1st 

Qk 

Du 

Jt 

it 

El 

it 

ik. 

El 

it 

it 

ikj 

it 

Ik. 

PM! 

it 

Jb 

m 

Ft.  Rots 

b9 

U9 

b9 

b9 

b9 

U9 

13 

35 

31 

b9 

b9 

b5 

b7 

U7 

U5 

Ft.  Collins 

E 

D 

b9 

U9 

U9 

b3 

b3 

1*3 

39 

1*5 

1*1 

b9 

U9 

U9 

37 

37 

37 

Dodg s  City 

27 

31 

25 

b9 

b9 

b9 

19 

23 

21 

27 

27 

27 

25 

25 

25 

1*5 

1*5 

U5 

Vicksburg 

31 

31 

31 

37 

35 

35 

U7 

U7 

1*7 

1*9 

1*9 

1*9 

b7 

b7 

b9 

31 

31 

31 

lin-flils 

15 

11 

13 

35 

35 

33 

b9 

U9 

1*9 

17 

19 

17 

27 

27 

25 

29 

1*3 

1*3 

Cairo 

15 

17 

11 

27 

25 

19 

15 

15 

27 

1*7 

1*7 

1*7 

27 

27 

27 

13 

13 

13 

Ihdison 

37 

b3 

b9 

31 

31 

33 

b7 

b7 

1*7 

U3 

1*5 

1*3 

29 

b7 

U7 

27 

27 

29 

Pittsburgh 

17 

U3 

U3 

39 

bl 

37 

9 

9 

1*9 

17 

bl 

bl 

19 

19 

19 

23 

b3 

23 

Iynchburg 

27 

27 

27 

U7 

b7 

37 

U5 

1*7 

1*7 

1*7 

1*5 

b7 

U9 

U9 

U9 

37 

7 

39 

STATION 

■I  | 

Pk 

,SBP 

siOkPlr 

2  OCT 

Sw  Qlr  Dir 

MOV 

St  Ok  Dir 

DEC 

d*.*. 

Ft.  Ross 

b9  U9  U5 

U  U  7 

37  37  b9 

31  31  b9 

b5  b5  b5 

bl  bl  bl 

Ft.  Collins 

U5  13  b7 

13  715 

17  17  33 

17  39  17 

39  b3  b5 

27  27  bl 

Dodg*  City 

1*3  U3  U3 

b5  1*5  1*5 

bl  bl  bl 

25  25  25 

39  bl  bl 

b5  b5  b7 

U3b3b3 

Vicksburg 

25  15  23 

17  21  21 

U3  1(3  It3 

b9  b9  b9 

17  17  13 

b3  b3  bl 

27  33  17 

Mtaphlf 

31  33  3$ 

b5b51*5 

b9  b9  1*9 

15  15  13 

b9  b9  U9 

33  33  33 

U7  b7b9 

Cairo 

23  27  31 

U  29  27 

27  27  27 

39  b9  b9 

U  9  21 

19  25  23 

17  25  23 

Ihdison 

U7  U7  U7 

1*3  39  39 

15  1513 

bl  bl  23 

13  U13 

b9  b9  b9 

b7  b7  b9 

Pittsburg* 

1*1  bl  U3 

23  23  25 

15  17  17 

19  27  5 

39  39  39 

19  19  23 

21  21  31 

313133 

lyncbburg 

b5  31  27 

31  33  3$ 

25  25  25 

19  19  23 

35  to  37 

5  5li7 
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biasing  this  sample,  as  in  Court's  study,  in  mean  and  in 
variance  as  depicted  on  the  bottom  of  Figures  5-10. 

The  main  difference  between  the  three  measures  is 
the  rate  of  change  as  k  becomes  larger.  shows  the 
greatest  amount  of  variance  while  shows  the  least. 

There  can  be  seen  no  significant  difference  between  and 
when  values  of  k  are  above  15.  With  biasing,  all  three 
behave  in  the  same  manner  and  a  detailed  discussion  of  the 
possible  significance  of  this  behavior  is  given  in  Court's 
paper. 

4.4  Analysis 

For  purposes  of  comparison,  the  seven  stations  for 
o 

which  S£m  values  have  already  been  determined  in  Court's 
study  (1968),  along  with  Fort  Collins  and  Fort  Ross  were 
chosen  for  study.  The  latter  two  stations  were  included 
so  that  data  from  drier  regions  might  also  be  studied.  Lo¬ 
cations  for  these  stations  along  with  the  number  of  years 
of  record  are  given  in  Map  20. 

In  a  preliminary  study  and  were  computed  for  the 
above  stations  and  the  values  of  k*  listed,  along  with  those 
for  S^,  in  Table  3.  Examination  of  this  Table  shows  that 
nc  method  of  computation  is  consistently  better  than  another 
and  that  in  the  majority  of  cases,  k*  values  are  identical. 
It  seems  obvious  that  such  an  approach  is  unsatisfactory 

for  the  purpose  of  determining  the  advantages  of  the  mean 
or  median. 
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AVERAGE  SQUARED  DIFFERENCE 
BETWEEN  THE  MEAN  FOR  K  YEARS 
AND  NEXT  YEAR'S  OBSERVATION. 


SAMPLE  3 

"•IOOO 

-.10 - .so— 


Figure  5.  (Nonna  1  Sample  No.  3,  Mated  in  mean,  Court,  1967,  Fig.  12) 


Figure  6.  (Normal  Sample  No.  3,  biased  in  variance,  Court,  1967,  Fig.  13) 
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AVERAGE  ABSOLUTE  DIFFERENCE 
BETWEEN  THE  MEAN  FOR  K  YEARS 
AND  NEXT  YEAR'S  OBSERVATION 
(BIASED  IN  MEAN). 


FIGURE  7. 
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YEARS o« RECORD 


A  second  table  was  then  prepared  (see  Table  4).  Here 

the  values  for  and  Dk  at  k*  are  given,  instead  of  the 

value  of  k*  itself.  As  mentioned  before,  is  not  com- 
2  2 

parable  with  S£  and  therefore  S£  is  not  included;  however, 
since  the  object  of  this  study  is  to  determine  the  relative 
merits  of  the  mean  and  median,  the  comparison  of  only 
and  Dfc  will  be  sufficient. 

Table  4  clearly  shows  that'  the  median  is  a  better 
value  than  the  mean  for  estimating  future  values.  In 
nearly  two-thirds  of  the  cases  studied  the  value  of  D^  is 
lower  than  that  for  . 

To  further  bring  out  the  advantages  of  the  median  over 
the  mean  Figures  11-13  were  prepared  showing  the  values  of 
Qk  and  Dk  for  the  odd  values  of  k  from  3  to  49.  It  can  be 
seen  that  a  fewer  number  of  years  is  required  to  obtain  Dk 
than  is  required  to  obtain  Qk  of  the  same  value. 

When  m  is  varied  from  one  to  ten  to  predict  beyond 
next  year,  it  is  found  that  k*  decreases  as  m  becomes  lar¬ 
ger.  This  phenomenon  occurs  regardless  of  what  statistic 
is  being  used.  This  "extension"  is  clearly  seen  in  Figure 
14.  The  reason  for  this  extension  is  not  clear  at  present, 
but  is  being  studied  by  Dr.  Court  and  will  be  included  in  a 
future  report. 
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TABLE  4,— OPTIMUM  VALUES  OF  Qk  AND  !\l 


STATION 

i  m 

‘J 

MAR 

Qlr  Dir 

APR 

Qlr  Dk 

MAI 

.  Sit.  ,9k 

mm 

Ft*  Ross 

9*96  8.21 

8.67  7.89 

6.1*9  5.81* 

U.88  1**89 

2 .91*  2.62 

1.78  1.1*7 

Ft*  Collins 

.61  .55 

•61*  .61 

1.15  1.15 

2.18  2.01 

2.79  2.71* 

2.12  2.17 

Dodgs  City 

•80  .71* 

1.06  1.07 

1.71  1.65 

2.82  2.73 

3.77  3.78 

3.70  3.57 

Vicksburg 

5.56  5.29 

5.09  5.09 

1*.38  U.38 

5.78  5.5U 

1*.1*8  1*.50 

l*.5l  l*.32 

Jfsaphis 

6.70  6.1*3 

1*.61*  1*.1*6 

U.99  U.97 

5.1*5  5.52 

1*.1*1*  1**65 

1*.15  U.15 

Cairo 

5.61*  5*51 

3.70  3.66 

5.^3  5.18 

3.1*3  3.56 

1**53  U.53 

5.13  U.91 

Mtdison 

1.87  1.86 

1.60  1.59 

1.71  1.72 

2.1*0  2.39 

3.55  3.1*1* 

3.33  3.33 

Pittsburgh 

2*72  2.67 

2.12  2.05 

2.81  2.83 

2.59  2.61* 

3.01  3.09 

3.33  3.39 

Iynchburg 

2.90  2.89 

2.36  2.29 

3.37  3.36 

2.30  2.30 

3.56  3.58 

3.89  3.91 

STATION 

JUL 

Qlr  Dir 

AUQ 

Qlr  Dk 

SEP 

Qlr  Dk 

OCT 

Qlr  Dir 

NOV 

Qlr  Dk 

1R 

Qlr  D|, 

Ft.  Ross 

.16  .13 

•27  .18 

1.58  1.12 

3.88  3J*1 

7.71  6.83 

9.71  9.19 

Ft.  Collins 

1.89  1.58 

2.35  2.18 

2.59  2.23 

2.05  1.96 

.78  .81 

.91  .81* 

Dodgs  City 

3.57  3.32 

2.75  2*68 

2.37  2.31* 

2.1*7  2.32 

1.59  1.60 

.91  .92 

11.88  11.81 

Vicksburg 

5.23  5.29 

3.96  3.99 

U.37  3.81* 

U.02  3.1*1 

5.31*  1*.S6 

1*.55  U.22 

18.11 18.13 

Memphis 

3.85  3.90 

3.52  3.25 

L.10  3.91 

L.38  U.29 

1*.1*2  U.27 

l*.8l  U.58 

22.81*  21.21 

Cairo 

3.78  3.93 

1*.32  1.32 

1*.02  1*.12 

3.68  3.66 

U.65  1*.72 

3.3U  3.36 

17.39  17.UU 

Msdison 

3.61  3.55 

3.62  3.59 

1*.50  U.39 

2.72  2.72 

2.1*8  2.51 

1.22  1.26 

9.68  9.76 

Pittsburgh 

3.33  3.16 

3.36  3.32 

3.06  2.97 

3.25  3.31 

2.53  2.1*8 

2.19  2.27 

10.30  10.30 

Iynchburg 

3.85  3.77 

5.21  5.01 

1*.75  U.33 

3.93  3.73 

3.15  3.06 

2.53  2.1*8 

15.12  15.07 

mi  mmn 

tor 
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AVERAGE  ABSOLUTE  DIFFERENCE  BETWEEN 
THE  MEDIAN  (A)  OR  MEAN  (B)  OF  K  YEARS 
AMD  NEXT  YEAR'S  OBSERVATION. 


DODGE  CITY 
Annual  Precipitation 


CHAPTER  V 


CONCLUSIONS 

In  climatology,  as  in  many  other  areas  of  geography, 
it  is  often  necessary  to  characterize  entire  sets  of  obser¬ 
vations  by  a  single  value.  In  the  past  the  arithmetic 
average,  or  mean,  has  been  the  most  generally  used  value. 
Many  variables  have  values  that  have  the  opportunity  for 
downward  variation  strictly  limited  by  zero,  while  upward 
variations  are  theoretically  unlimited.  Yet  it  is  the 
lower  values  that  dominate  the  geography  of  an  area.  Here 
the  use  of  the  mean  is  particularly  at  fault. 

The  arithmetic  mean  has  been  shown  by  the  many  maps 
contained  in  this  thesis  to  be  a  questionable,  or  at  least 
misleading,  descriptive  value.  The  arithmetic  mean  has 
been  hitherto  employed  as  a  measure  of  central  tendency  for 
convenience  of  description.  The  indubitable  advantages  of 
the  arithmetic  mean  are  the  ease  with  which  it  is  calcula¬ 
ted  (at  least  before  the  advent  of  the  computer),  and  the 
fact  that  the  sum  of  the  twelve  monthly  normals  is  equal  to 
the  average  of  annual  totals.  Its  chief  fault  lie"  in  the 
fact  that  equivalence  of  weight  is  given  to  each  unit  by 
which  an  individual  record  departs  from  the  mean. 

Climates  ar<=  not  c.onstant  in  position  and  it  ir  falla¬ 
cious  to  regard  them  as  being  so.  The  search  for  any  valid 
mean  expression  would  seem  foredoomed  to  failure.  Climatic 
statistics  must  therefore  be  examined  in  their  entirety, 
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and  obsolete  normals  must  be  replaced  by  indices  of  proba¬ 
bility,  as  exemplified  in  the  Atmospheric  Humidity  Atlas-- 
Northern  Hemisphere  (Gringorten  et  al_,  1966)  .  The  first 
step  toward  this  goal  would  be  to  replace  the  arithmetic 
mean  with  the  median. 

Many  studies  are  liable  to  the  biased  errors  inherent 
in  the  arithmetic  mean  and  the  temptation  to  beg  the  whole 
question  when  phenomena  are  submitted  to  subjective  examin¬ 
ation.  Neither  of  these  dangers  would  seem  to  offer  much 
difficulty  when  squarely  faced;  but  the  more  elaborate  and 
refined  the  subsequent  analysis,  the  easier  it  is  to  over¬ 
look  initial  limitations. 

The  present  notion  that  the  mean  value  is  something 
that  will  recur,  or  that  it  is  the  value  which  best  repre¬ 
sents  what  we  expect  to  happen  again,  should  be  replaced. 
The  knowledge  that  the  mean  is  only  one  measure  of  central 
tendency,  and  a  poor  one  at  that  when  considering  many  of 
the  climatological  variables,  supports  such  replacement .  In 
the  future  it  is  hoped  that  summaries  of  climatological  data 
will  feature  the  median  along  with,  or  instead  of,  the  mean. 

Climatology  is  not  the  only  subfield  of  geography 
where  skewed  distributions  are  encountered.  Many  geographi¬ 
cal  distributions  are  singly  or  doubly  bounded  and  since 
geographers  are  interested  in  reducing  large  amounts  of  data 
to  representative  values  which  can  then  be  used  to  describe 

an  area  or  compared  with  other  values,  it  is  of  the  utmost 
importance  that  the  values  are  comparable  and  best  estimate 


9 


the  actual  distribution.  While  the  moan  may  have  some 
computational  advantages,  it  still  fails  to  represent  the 
most  typical  value  in  many  distributions.  This  failure 
has  been  recognized  for  a  few  extremely  skewed  distribu¬ 
tions,  and  the  median  values  are  now  in  general  use  to 
represent  ’average'  income  and  'average'  number  of  years 
of  education.  But  the  arithmetic  mean  < ontinues  to  be  used 
to  portray  values  of  agricultural,  indur trial,  and  mineral 
production;  to  indii.  .  te  consumption  of  goodr;  to  summarize 
climatic  conditions;  and  many  other  variables  which  might 
be  better  portrayed  by  the  median. 
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APPENDIX  I 


TABLE  1. —PRECIPITATION  STATIONS  FOR  WHICH  MONTHLY  VALUES  ’WERE  USED 

(Discussed  on  Pages  10*13) 


SEQUENCE 

NUMBER 

STATION 
I.D.  NO. 

STATION  NAME 

LATITUDE 

NORTH 

LONGITUDE 

WEST 

YEARS  OF  RECORD 
(ENDING  I960) 

1 

2.185 

Clifton 

33.05 

109.28 

52 

2 

2.359 

Grand  Canyon  National  Park 

36.05 

112.12 

1*0 

»  3 

2.571* 

Mount  Trumball 

36.  h2 

113.33 

3C 

1* 

2.632 

Pays  on  RS 

3U.23 

111.33 

31* 

5 

2.656 

Pinal  Ranch 

33.35 

110.98 

65 

6 

2.680 

Prescott 

31*.  55 

112.1*5 

52 

7 

2.882 

Tucson  University  of  Arizona 

32.23 

110.95 

66 

8 

2.965 

Yuma  Citrus  Station 

32.62 

111*  .65 

1*0 

9 

3.023 

Arkansas  City 

33.62 

91.20 

72 

10 

3.01,6 

Bates ville  Land  D.  No.  1 

35.75 

91.63 

61 

11 

3.160 

Conway 

35.08 

92  .*,7 

77 

!  12 

3.21*1* 

Fayettesville  Exp.  Station 

36.10 

91*.  17 

70 

13 

?.1»76 

Mena 

3l*.58 

9i*.25 

50 

11* 

3.5ol* 

Mountain  Home  1  NNW 

36.33 

92.38 

uh 

15 

3,582 

Pochontas 

36.27 

90.98 

67 

16 

3.693 

Subiaco 

35.30 

93.65 

63 

, 

U.023 

Antioch  F.  Mills 

38.02 

121.77 

31 

'  18 

1*.038 

Auburn 

38.90 

121.07 

61 

19 

L.076 

Big  Creek  Power  House 

37.20 

119.25 

1*5 

'  20 

U.079 

Big  Sur  State  Park 

36.25 

121.78 

1*6 

fcl 

U.170 

Chester 

1*0.30 

121.22 

50 

22 

U.316 

Fort  Bragg 

39.95 

123.30 

61 

23 

U.319 

Fort  Ross 

38.52 

123.25 

85 

21* 

1*.1*02 

Hollister 

36.85 

121.1*0 

87 

25 

i  ■. 

U.522 

Lytle  Creek  Power  House 

3U.20 

117.1*5 

55 

26 

1*.51*5 

McCloud 

1*1.27 

122.13 

50 

27 

U.612 

Needles 

31*.  77 

111*. 62 

69 

28 

L.618 

Newport  Beach  Harbor 

33.6C 

117.88 

30 

29 

U,.  61*0 

OJai 

31*  .1*5 

119.25 

56 

30 

U*.  771* 

San  Diego  WB  APT. 

32.75 

117.17 

111 

75 


31 

32 

33 

31* 

35 

1*.78$ 

l*.80l* 

U*835 

1*.897 

1**9014 

San  Lula  Obispo  Poly 

Scotia 

Sonora 

Topanga  Patrol  Station  FC  6 
Trona 

35.30 

j*0.1*8 

37.98 

3l*.08 

35.78 

120.67 

121*.10 

120.38 

118.60 

117.38 

91 

35 

73 

30 

1*1 

36 

37 

38 

39 

1*0 

1**909 

1*.913 

1*.91*5 

i*.91*9 

1**970 

Tustin  Irvine  Ranch 

Twin  Lakes 

Wasco 

Weavervllle  Rs. 

Willows 

33.73 

38.70 

35.60 

1*0.73 

39.53 

117.78 

120.05 

119.33 

122.93 

122.20 

81* 

38 

61 

71 

82 

la 

1*2 

1*3 

1*1* 

1*5 

5.129 

5.153 

5.156 

5.218 

5.21*3 

Cannon  City 

Cheesman 

C!  jyenne  Wells 

Del  Norte 

Durango 

38.1*3 

39.22 

38.82 

37.67 

37.26 

105.27 

105.28 

102.35 

106.35 
107.88 

67 

58 

61* 

31* 

66 

1*6 

1*7 

1*8 

1*9 

50 

5.300 

5.301* 

5.1*l*l 

5.1*83 

5.572 

Fort  Collins 

Fort  Morgan 

Julesburg 

Las  Animas 

Montrose  No.  2 

1*0.58 
1*0.25 
ia.  00 
38.07 
38.1*8 

105.08 

103.80 

102.25 

103.22 

107.88 

63 

51* 

1*9 

91* 

70 

51 

52 

53 

51* 

55 

5.762 

5.791* 

5.929 

10.001 

10.01*5 

Shoshone 

Steamboat  Springs 

Yuma 

Aberdeen  Expt.  Station 
Arrowrock  Dam 

39.57 

1*0.50 

1*0.12 

1*2.95 

1*3.60 

107.23 

106.83 
102.73 

112.83 
115.92 

51 

51 

71 

1*6 

1*9 

56 

57 

58 

59 

60 

10.11*1 

10.271 

10.391* 

10.501 

10.651* 

Cambridge 

Dubois  Expt.  Station 

Hailey  RS 

Kooskia 

Oakley 

1*1*. 5  7 
1*1*.25 
1*3.52 
1*6.15 
1*2.23 

116.68 

112.20 

111*  .32 
115.98 
113.88 

58 

39 

52 

52 

67 

61 

62 

63 

61* 

65 

10.808 

10.811* 

13.036 

13.221 

13.523 

Salmon 

Sandpoint  Expt.  Station 
Atlantic  1  NE 

Des  Moines  WB  City 

Mason  City  3N 

1*5.18 
1*8.28 
ia.i*  2 
1*1.58 
1*3.16 

113.88 

116.57 

95.00 

93.62 

93.20 

38 

50 

70 

83 

58 

66 

67 

68 

69 

70 

13.639 

13.716 

11*. 177 
U*.187 

11*  .21*6 

Ottuma 

Rockwell  City 

Concordia  WB  City 

Council  Orove 

Ellsworth 

1*1.00 

1*2 .1*0 
39.57 
38.67 
38.73 

92  .1*3 

91**62 

97.67 

96.50 

98.23 

68 

65 

75 

52 

56 

76 


71 

15.376 

Holton 

72 

15.552 

La  Cygne 

73 

15.517 

Medicine  I-ogrie 

71* 

15.637 

Phillipsburg 

75 

15.653 

Plains 

76 

1II.66U 

Quin ter 

77 

15.731 

Sedan 

78 

15.731 

Sedgwick 

79 

15.819 

Toronto 

80 

16.11a 

Calhoun  Expt.  Station 

81 

16.570 

Jenningr 

82 

16.612 

Melville 

83 

16.666 

New  Orleans  WB  City 

85 

16.735 

Plain  Dealing 

85 

16.892 

Tallulah  Delta  Lab. 

86 

23.130 

Capringer  Mills 

87 

23.158 

Chillicothe  25 

88 

23.223 

Dexter 

89 

23.250 

Eldon 

90 

23.282 

Fayette 

91 

23.3  Oli 

Frederick town 

92 

23.379 

Hermann 

93 

23.598 

Neosho 

95 

23.772 

Shelbina 

95 

23.871 

Warrensburg 

96 

23.899 

Willow  Springs 

97 

25.036 

Augusta 

98 

2U.CU3 

Sallnntine 

99 

25.077 

Pig  Sandy 

100 

25.105 

Bozeman  Agri,  College 

101 

25.260 

East  Anaconda 

102 

25.269 

Ekalaka 

103 

2U.31U 

Fortine  Inne 

loa 

25.389 

Hamilton 

105 

25.398 

Haugan 

106 

25.552 

Jordan 

107 

25.529 

Lustre  5  NNW 

108 

25.576 

Moccasin  Expt.  Station 

109 

25.729 

Saint  Ignatius 

no 

25.881 

West  S lacier 

39.1*7 

95.73 

58 

38.35 

95.77 

30 

37.27 

98.58 

68 

39.77 

99.32 

69 

37.27 

100.58 

51 

39.07 

100.23 

30 

37.12 

96.17 

76 

37.92 

97.53 

55 

37.8C 

95.95 

65 

32.52 

92.33 

69 

30.23 

92.67 

63 

30.68 

91.75 

75 

29.95 

90.07 

91 

32. 90 

93.68 

67 

32.  ho 

91.22 

53 

37.80 

93.80 

35 

39.75 

93.55 

53 

36.80 

89.97 

37 

38.35 

92.58 

57 

39.15 

92.68 

76 

37.57 

90.30 

36 

38.70 

91.55 

96 

36. 87 

95.37 

73 

39.68 

92.05 

•  7 

38.77 

93.73 

77 

36.98 

91.97 

37 

57.58 

112.38 

35 

1*5.95 

108.13 

51 

59.17 

110.12 

37 

55.67 

113.05 

58 

56.10 

112.92 

55 

5?. 97 

l0u.;-'3 

56 

58.78 

115.90 

39 

■'*6.25 

115.15 

33 

57.38 

115.50 

58 

57.32 

106.90 

30 

58,55 

105.93 

39 

57.05 

109.95 

30 

57.32 

115.10 

52 

5  8.5c 

113.99 

35 

1 


111 

25.093 

Blair 

1*1.55 

96.13 

91 

112 

25.115 

Bridge  Port 

1*1,67 

103.10 

63 

113 

25.202 

Crete 

1*0,62 

96.95 

81 

11 It 

25.280 

Ewing 

1*2.25 

98.35 

68 

115 

25.302 

Port  Robinson 

1*2.67 

103.1*7 

77 

116 

25.318 

Genoa 

1*1.1*5 

97.73 

85 

117 

25.363 

Hartington 

1*2.62 

97.27 

69 

118 

25.697 

Purdum 

1*2.07 

100.25 

58 

119 

25.701* 

Ravenna 

1*1.03 

98.92 

83 

120 

26.005 

Adaven 

38.12 

115.58 

1*2 

121 

26.257 

Elko  >13  Apt. 

1*0.83 

115.78 

91 

122 

26.517 

Ming 

38.38 

118.10 

53 

123 

26.678 

Reno  WB  Apt. 

39.50 

119.78 

90 

12U 

29.152 

Carrisoso 

33.65 

105.88 

92 

125 

29.181 

Cimarron 

36.52 

Id*.  92 

57 

126 

29.191* 

Clovis 

31*  .1*0 

103.20 

1*9 

127 

29.285 

Elida 

33.95 

103.65 

1*5 

128 

29.327 

Fort  Bayard 

32.80 

108.15 

90 

129 

29.378 

Hachita 

31.92 

108.32 

51 

130 

29.1*71* 

lake  Avalon 

32.U8 

10l*.25 

1*6 

131 

29.666 

Pecos  RS 

35.58 

105.66 

37 

132 

29.851* 

State  University 

32.28 

106.75 

100 

133 

29.990 

Zuni  FAA  AP. 

35.10 

108.78 

1*6 

13b 

32.219 

Dickenson  Expt.  Station 

1*6.88 

ioe.8o 

69 

135 

32.362 

(brand  Porks  U. 

1*7.92 

97.08 

69 

136 

32.1*1*2 

Jamestown  St*  Hosp* 

1*6.86 

98.68 

66 

137 

32.561* 

Max 

1*7.82 

101.30 

30 

138 

32.602 

Mohall 

1*8.77 

101.52 

67 

139 

3U.350 

Oeary 

35.63 

98.32 

1*9 

UiO 

31*.1*1*5 

Idabel 

33.90 

91*. 82 

31* 

11*1 

3l*.l*77 

Kenton 

36.92 

102.97 

60 

11*2 

31*. 693 

Pauls  Valley 

3U.75 

97.22 

57 

11*3 

31*.  701 

Perry 

36.28 

97.28 

30 

11*1* 

31*. 91*5 

Webber  Falls 

35.52 

95.13 

62 

1 1*5 

31*. 963 

Wichita  Mt.  Wlr 

31*. 73 

98.72 

1*5 

11*6 

35.020 

Antelope  IN 

1*1*.  92 

120.72 

36 

11*7 

35.069 

Bend 

1*1*.07 

121.32 

50 

11*8 

35.190 

Cottage  Or ore  1  S 

1*3.78 

123.07 

1*1* 

11*9 

35.211* 

Danner 

1*2.93 

117.33 

30 

190 

35.269 

Estachada  2  SE 

1*5.27 

122.32 

52 

78 


r 


151 

152 

153 
152* 
155 

35.3UU 

35.383 

35.1*67 

35.561 

35.691 

Grants  Pass 

Heppner 

Lakeviev 

Minam  7  NE 

Prospect  2  SW 

1*2. 1*3 
1*5.33 
1*2.  ie 
1*5.68 
1*2.73 

123.32 

119.55 

120.35 

117.60 

122.52 

72 

55 

1*8 

51 

53 

156 

157 

158 

159 

160 

35.725 

35.905 

39.030 

39.197 

39.280 

Rock  Creek 

Warn  Springs  Reservoir 
Amour 

Cotton  Wood 

Eureka 

1*1*.75 

1*3.57 

1*3.32 

1*3.97 

1*5.77 

118.08 

118.20 

98.35 

101.87 

99.62 

Ul 

30 

63 

51 

52 

161 

162 

163 

16U 

165 

39.383 

39.U01 

39.1*66 

39.1*86 

39.551* 

Highmore  1  W 

Hot  Springs 

Ladelle  7  NE 

Lennon 

Milbank 

1*1*.  52 
1*3. 1*3 
1*1*  .68 
1*5.93 
1*5.22 

99.1*7 

103.1*7 

98.00 

102.17 

96.63 

58 

57 

6U 

Ul* 

71 

166 

167 

166 

169 

170 

39.767 

39.855 

39.91*1* 

1*1.012 

1*1.050 

Sioux  Fall  WB  AP 

Vale 

Wood 

Albany 

Balmorhea  Exp.  Station 

1*3.57 
1*1*  .62 
1*3.50 
32.73 
31.00 

96.73 

103.1*0 

100.1*8 

99.30 

103.68 

70 

52 

U8 

79 

37 

171 

172 

173 

17i* 

175 

1*1.061 

1*1.111* 

1*1.202 

1*1.318 

1*1.31*3 

Beaumont 

Brownwood 

Corsicana 

Flatonia 

Galveston  WB  City 

30.08 

31.72 

32.03 

29.68 

29.30 

91*. 10 

98.98 

96.1*7 

97.10 

91*.  83 

68 

68 

75 

53 

89 

176 

177 

178 

379 

3  80 

1*1.351 

1*1.373 

1*1.1*08 

1*1.1*78 

1*1.50? 

George  West 

Greenville  2  SW 

Henderson 

Kerrville 

Lampasas 

28.35 

33.12 

32.15 

30.03 

31.05 

98.12 

96.13 

91*.  80 

99.13 

98.18 

L5 

59 

52 

65 

66 

181 

182 

183 

18U 

185 

1*1.597 

1*1.695 

1*1.721 

1*1.726 

1*1.765 

Mission 

Perry  ton 

Post 

Presidio 

Riverside 

26.22 
36.1*0 
33. 2C 
29.55 
30.85 

98.32 

100.  G2 

101.37 

101*. 1*0 

95.1*0 

LO 

U8 

UU 

30 

57 

186 

187 

188 

189 

190 

1*1.863 

1*1.928 

Ul. 933 
1*1.953 
1*2.210 

Sterling  City 

Valley  Junction 

Vega 

Weatherford 

Deseret 

31.85 

30.83 

35.25 

32.75 

39.28 

100.98 

96.63 

102.1*3 

97.80 

112.65 

30 

58 

30 

67 

61 

79 


u 


191 

1*2.300 

Port  Duchtsne 

1*0.28 

109.85 

71 

192 

1*2.390 

Hiawatha 

39.1*8 

111.02 

39 

193 

U2.U51 

Kanab  Power  Houtt 

39.05 

112.52 

39 

19U 

U2.515 

Loa 

38.1*0 

111.65 

57 

195 

1*2.565 

Milford  MB  Apt. 

38.1*3 

113.02 

1*1* 

196 

1*2.727 

Richmond 

1*1.90 

111.82 

1*9 

197 

1*2.812 

Spanish  Fork  Power  Hooss 

1*0.08 

111.60 

50 

198 

1*2.877 

Tooala 

1*0.53 

112.30 

61* 

199 

1*5.092 

Brooklyn 

1*6.77 

123.52 

30 

200 

1*5.122 

Cedar  Lake 

U7.1*2 

121.95 

58 

201 

1*5.135 

Chelan 

1*7.83 

120.03 

69 

202 

1*5.159 

Colfax  1  Mf 

1*6.88 

117.38 

69 

203 

1*5.322 

Qoldendale 

1*5.82 

120.83 

50 

201* 

1*5.355 

Hatton  8  E 

1*6.77 

118.67 

56 

205 

1*5.1*77 

Longview 

1*6.17 

122.92 

36 

206 

1*5.581* 

kewhalem 

1*8.68 

121.25 

36 

207 

1*5.701* 

Rimrock  Teton  Dan 

1*6.65 

121.13 

51 

208 

1*5.751 

Sedro  Wolley  1  E 

1*8.50 

122.22 

61* 

209 

1*5.758 

Shelton 

1*7.20 

123.10 

30 

210 

1*5.821 

Sunnyaide 

1*6.32 

120.00 

66 

7T 

211 

1*5.833 

Tatoosh  Island  WB 

1*8.38 

12U.73 

77 

212 

1*5.938 

Winthrop  1  WSW 

1*8 .1*7 

120.18 

39 

213 

1*8.118 

Buffalo  Bill  Dam 

1*1*. 50 

109.18 

1*8 

211* 

1*8.272 

Dubois 

1*3.55 

109.62 

1*8 

215 

1*8.1*07 

Green  River 

1*1.53 

109.1*8 

1*0 

216 

1*3.583 

Lusk 

1*2.77 

10U.U3 

62 

217 

1*8.710 

Pathfinder  Dam 

1*2 .1*7 

106.83 

60 

218 

1*8.816 

Sheridan  Field  Station 

1*1*.85 

106.87 

1*3 

219 

1*8.991 

Yellowstone  Park 

1*1*.97 

110.70 

72 
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CLIMATIC  PREDICTION 


Optimum  Length  of  Record 

1#  Origin  of  Program 

This  program  was  prepared  in  the  Climatology  Laboratory  of  San  Fer¬ 
nando  Valley  State  College,  Northridge,  California  by  Paul  E.  Roy,  Jr.  and 
William  F.  Slusser  as  part  of  Air  Force  Project  "Optimum  Record  for  Clima¬ 
tic  Estimation  and  Prediction,"  to  assist  in  the  analysis  of  climatic  data 
to  determine  the  optimum  length  of  record  for  climatic  prediction. 

2.  Purpose  of  Program  > 

This  program  accepts  monthly  and  yearly  data  in  varied  formats  (Sec.  5) 
and  computes  the  extrapolation  variance,  S(k,m),  the  absolute  prediction 
error  using  the  mean,  Q(k,ra),  and  the  median,  D(k,m,*  for  varying  k  year 
periods  and  observations  m  years  ahead.  It  then  finds  the  optimum  length 
of  record  (k*)  for  each  of  these  and  computes  F(k,ro),  QF(k,ro),  and  DF(k,m), 
the  percentage  difference  between  the  value  obtained  at  k*  and  that  found 
using  the  other  values  of  k. 

3.  Description  of  Equipment 

This  program  was  developed  for  use  on  the  IBM  360/75*  but  can  be  used 
on  most  computers  using  Fortran  IV* 

Method  of  Computation 

a*  S(k,ra)  is  obtained  by  averaging  the  squared  differences  between 
the  mean  of  k  successive  observations  and  an  observation  ra  ye?.rs  later  for 
which  the  value  is  to  be  estimated  from  that  for  the  k  years. 
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Q(k,ra)  and  D(k,m)  are  obtained  by  averaging  the  absolute  differences 
of  the  same  periods* 

Values  of  k  from  1  to  any  desired  length  of  record  nay  be  used  to 
obtain  values  for  S(k,m)  and  Q(k,m)  while  only  odd  values  of  k  beginning 
with  3  0^0  used  for  finding  values  of  D(k,m). 

Lastly ,  a  varies  from  1  to  10* 

b.  F(k,m),  QF(k,m),  and  DF(k,m)  are  obtained  by  dividing  all  values 
of  S(k,m),  QF(k,m),  and  D(E,m)  respectively  by  the  minimum  vain  obtained 
for  each  and  multiplying  by  100*  These  minimum  or  optimum  values  are  Indi¬ 
cated  on  the  print  out  by  an  asterisk  (*),  and  the  value  of  k  at  which  the 
value  is  obtained  is  known  as  the  optimum  length  of  record  (k»)* 

Sums  for  F(k,l),  QF(k,l),  and  DF(k,l)  of  the  twelve  monthly  values 
for  each  value  of  k  are  given  on  the  last  page  of  output. 

c.  The  program  also  supplies  a  print  out  of  all  data  used  in  the 
calculations  serving  as  a  check  to  see  that  the  data  is  properly  read. 

5*  Input 

The  format  of  the  data  to  be  examined  is  read  In  so  that  the  program 
is  readily  usoable  (see  Input  description).  The  program  uses  monthly  and 
annual  data  in  chronological  order.  Since  a  doubly  dimensioned  array  is 
used,  each  card  should  contain  all  or  any  number  of  months  of  one  year  only, 
along  with  the  annual  value  if  desired.  The  order  of  the  months  on  the  data 
card  must  coincide  with  the  order  of  the  months  on  card  B  (see  Input  descrip¬ 
tion)* 

6.  Output 

All  output  is  in  printed  form,  as  summarized  by  the  illustration* 
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7»  Operating  Instructions 

Standard  Fortran  IV  is  us ad.  No  sense  switches  are  used.  At  the  pre¬ 
sent  the  program  does  not  handle  tapes,  but  it  can  be  easily  modified  .to  do 


8»  Description  of  Terms 

IN  -  physical  input  unit 

IOUT  -  physical  output  unit 

NST  •  the  symbol  for  denoting  the  minimum  value 


XWON(I) 


FMT(I) 

X(L,NH) 

SIIF(K) 

SAF(K) 

'SF(K) 

XMIN(M) 

XMIN(M) 

KJ) 


SSUfflL 


SSUM 


S(K,M) 


-  the  symbol  denoting  other  than  minimum  value 

-  thirteen  three  letter  words:  Jan,  Feb,  •  •  .Sum;  or  Ir 


-  the  maximum  value  of  k 


N  -  the  number  of  observations  per  month  (or  year) 

ISETS  -  the  number  of  sets  of  data  por  card 

NNAME(l)  -  the  name  of  the  station,  the  variable,  and  the  years 


of  record 


-  the  format  of  the  data  being  read  in 

-  value  of  observation,  monthly  or  annual 

-  sum  of  monthly'  values  of  DF(k,m) 

-  sum  of  monthly  values  of  QF(k,m) 

-  sum  of  monthly  values  of  F(k,m) 

-  minimum  value  of  Q(k,m) 

-  minimum  va  lue  of  S(k,m),  SQ  or  D(k,m) 

-  the  mean  or  median  value  of  a  k  year  period 

-  sum  of  absolute  differences  using  mean 

-  sum  of  absolute  differences  using  median  or  squared 
differences  ii^ng  mean 

-  extrapolation  variance  or  absolute  prediction  error 


using  mean 
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Q(K£M)  •  the  absolute  prediction  error  using  the  man 

LOC(M)  -  the  value  of  k  where  NST  la  to  be  placed 

NBIANK(K,M)  •  the  symbol  to  be  printed  after  each  value  of  S(k,m), 

SQ  or  D(k,a) 

MBIANK(K,M)  -  the  symbol  to  be  printed  after  each  value  of  Q(k,m) 

Z(MID)  -  the  median  of  k  Tears 


INPUT  DESCRIPTION 

A*  One  card  containing  an  asterisk  (*)  in  the  first  column* 

B*  One  card  containing  the  names  of  the  months  using  columns 
1  through  36,  three  columns  per  month*  If  annual  data  is 
Included,  sum  or  yr  should  be  punched  in  columns  37  through 

39. 

C*  One  card  containing  the  maximum  value  of  k  desired  in 
columns  1  through  3 I  the  number  of  years  of  record  in 
columns  U  through  6;  the  number  of  sets  of  observations  (13 
if  all  months  and  the  annual  amount  are  included)  in  columns 
7  through  8;  and  the  station  name,  variable  name,  and  date 
of  beginning  and  end  of  observations  in  the  remainder  of  the 
columns* 

D*  One  card  containing  the  format  of  the  data  in  the  first  60 
columns* 

2*  Data 

F*  Each  successive  group  of  data  requires  only  cards  C  and  D* 
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COMPUTER  PROGRAM  ORCEP  (FORTRAN  IV,  IBM  360/75)  Pfcgd  1 

_ DIMENSION  X<  13,206),  Y(  206  ).SF  (50)  ,S(50  , 1  1  )_,J_(  5J  )  ,  NBI.  ANK(50, 11  )  , 

1 XM I N ( 11  ),LOC( 11),XM0N( 1 3 ) , NT  BL ( 50 ) ,  NNAME  (  24  ) ,  SUM  ( 1C~0 ) 

_ D1MLNS1CIN  Z(  100)  ,SAF(5C)  ,Q(50,  ID,  M  BLANK  (50,11),  LOCI  (_ID _ 

DIMENSION  FMT(20),XMiNl(ll),SDF(50) 

_ LfORMAT  ( //// 1  AH f(K,M)  FOP  ,A3,5H  OF  A,14,13H  YEAR  PERIOD  ,24A3 

1 ,  /// 3H  M-  ,  6 X  » lFl ,  3X,  1H2 , 8X,  1H3,  «X  ,  1H4,  8X ,  i  H5 , 8X ,  1H6, 8X,  1H7,  8X,  1H8, 
2BX, 1H9,  7X,2H10, _ /3H  K,  /  ) _ 

2  FORMAT  (////14H  F  (  K »  M  )  FOR  ,A3,5H  OF  A, 14, 13H  YEAR  PERIOO  ,24A3 

_ l,_///7X,4HMEAN,JX,13HABS0LUTc  ME  AN ,  IX  ,6HMED  I  AN, /3H  K,/) 

3  FORMAT  <////////  10H  DATA  FOR  ,A3,5H  OF  A,  14,13H  YEAR  PERIOD  ., 

. _ 124A3,//,21(  1H  ,1  OF  10.1,//)) 

4  FORMAT  ( / / / / 1  4  H  S(K,M)SQ  FOR  , A3 , 5H  OF  A,I4,13H  YEAR  PERIOD  ,24A3 

I, ///3H  M=,6X,  IM1,8X, 1H2,8X,1H3,8X, 1 H4 , 8X , 1 H5 , 8X , 1H6,8X, 1H7,8X, IH8, 
28X,1H9,7X, 2H10,  /3H  K,/) 

_  5  FORMAT  (  1H  ,J3 , 1  X,  1 2  ( F  8. 2,  A1  )i 
6  FORMAT  (1H  ,I3,ix,3F8.3) 

_ 1  FORMAT  (13A3J _ 

8  FORMAT •  02A1  ) 

9  FORMAT  (  IH1  ) 

10  FORMAT  <////15H  QF ( K , M )  FOR  ,A3,5H  OF  A , I 4 , 1 3M  YfAR  PERIOD  ,24A3 

_  1, ///3H  M=,6X,1H1 ,8X, IH?,8X,IH3,8X, IH4,8X,1H5,CX,  1 H6 , 8X , 1H7 , 8X ,  1H8, 

28X , IH9,  7X , 2H 1 C ,  /3H  K , / ) 

11  FORMAT  ( /// / 1 5H  Q(K,M)  FOR  ,A3,5H  PF  A, 14, 13H  YEAR  PERIOD  ,24A3 

1, ///3H  M=,6X, IH1,8X, IH2 , 8X , 1 H3 , 8X  ,  1H4 , 8X , 1 H5 , 8X , IH6, 3X , 1H7 , 8X, 1H8 , 
28X,IH9,7X,?H10, _ /3H  K,/J  _ _ 

12  FORMAT  (  /// / 1 5H  DF(K,M)  FOR  , A3 , 5H  OF  A, 14, 13H  YEAR  PERIOD  j?4A3 

1,  ///3H  Ms,6X,j|Hl_,3X,  1H2,8X,1H3,8X,  1H4,8X,  1H5,8X,  1H6,  SX,  IH7.8X,  1H8* 
28X,1H9,7X,2H10,  /3H  K,/) 

14  FORMAT  ( /// / 1 4H  D(K,M)  FUR  ,A3,5H  OF  A,I4,13H  YEAR  PERIOD  ,24A3 
1, ///3H  M=,6X,  1H1 ,8X, 1H2,8X,IH3,8X,  IH4,8X,1H5,8X, 1H6,RX, 1H7,8X,1H8, 

_ 28X,1H9,7X,2H10, _ /3H  K*/J _ _ 

15  FORMAT  (213, 12.24A3) 

16  FORMAT  (20 A3) _ 

I  N  =  5 

I 0UT  =  6  _ 

RE  AD  ( IN,  8  )N'S  T  ,  NBL 

_ READ(1N,7)(XM0N(I),1  =  1,13) _ 

DG  53  I XK= 1 , 7 

C  READ_IN  DATA _ 

READ( IN, 15 )KMAX,N, ISETS.NNAME 

READJ  IN,  16)  FMT _ 

DO  17  N  H  =  1 » N 

_  RE  AO ( IN , FMT ) ( X (L  ,NH ) , L  =  D I SE  T  S) 

17  CONTINUE 

DC  18  K  =  DKMAX  _ ____ _ _ 

SDF ( K)  =  0. 0 

SAF(K)  =  0._0  _ _ 

'  18  SF ( K  )  =  6  • 0 

C  CALCULATES  DEPARTURES  FROM  MEAN 
DO  51  I  1  =  1, ISETS 
LL=  1 1 
X1=II 

DO  19  M= 1,11  _ 

XM  INI  (M) =9999999. 

_ 19  XX  1  N ( M )  =  9999999 . _ 

WRi  IE( I  OUT, 9) 

_■  WRITf  (IC'UT,  3)XMGN(  D  )  ,  N  ,  NNAMf  ,  (  X  ( 1.  L  ,  NH )  ,  NM=  1 ,  N I 
. WRITL(inUT.9)  "  85 
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COMPUTER  PROGRAM  ORCEP  (FORTRA*  IV.  IBM  360/75) 


_ El  1).-0«Q 

DO  27  K=l,KMAX 

_ L=N-K _ 

C5=L*l 

_ C2=1.ZC5 _ 

C4=K 

_ Q1.ZC4 _ 

C 1=  1  •/ !C4*C4) 

_ m  l=Yl iitXtLLiK) _ 

DO  20  MM=1 «  L 

_ _ 

NN=MM+K 

_ Y ! J J I =Y (  MM  I+!X(LL«  NN  )-X!LL«MMl) 

20  CONTINUE 

_ D0._26K=LilQ _ _ 

SSUML=0.0 

_ SiLlMgjQaQ _ : _ 

NMK  I  =  N-M-K-f  1 

_ DO  21  J=1 tNMKl 

JKMl3  J+K+M-l' 

_ S.SUM  L=S  SUMl+LY!  JI *C-XI LL. JKM 1 1 1 ♦( Y ! J»*C-X(LL. JKM1 II 

21  $SUM=SSUM+ABS( Y( J) *C-X( LL» JKM1 ) ) 

_ C.6=NMKJL _ 

S!K,M)=SSUM1/C6 

_ Q  IK«M  )sSSlJM/Cfe _ _ _ 

IF !S!K,M)-XMIN1! MI  122,23,23 

22—XM1N  l  IMJ  =  S  !  K ,  M  _ 

LOCI I M  > =K 

23_MBLANK!K,M)=NBL  _ 

IF ( Q( K,M)-XMIN ( M) ) 24,25,25 

24  XM  IN(Hl=Q  (K # Ml _ _ 

LOC (M)=K 

25.  NBLANK.iK,M)  =NBL  _ 

26  CONTINUE 

27_C0NIINUE _ 

DO  28  M»l,10 

_ KKl=lflCltMl _ 

KK=LOC(M) 

_ MBLANK!KK1,M)=NST _  _ 

28  NBLANK! KK,M ) =NST 

WRITE!IQUT,4)XP0N( II )  ,N,NNAME..  .......... _ 

DO  29  K= 1 »  KMAX 

29  HRITE(IQUI-«.5J.K.»l.S.l  K, M) tMBLANKIKtMI  «Mrl.«lQl 
WRITE! IOUT , 9 ) 

_ HP  Lit  ( IOUT  ,-LLlXMON.LLI  1_*N,  NNAME _ 

DO  30  K= l , KMAX 

30L.WR  1.TE  !  I QUI ,5JK ,  IQiJLtNJ  iN  BLANK  IKaM). tK=_l  xlQ  i 
WRITE! IOUT , 9 ) 

HR  I T EC  LOUT, 1 IXMON! II) «N«NNAME _ 

00  33  K=  1 , KMAX 

_ DQ__3 1  J!=l  tlQ _ 

S(K,M)=<S!K,M)/XM INI (MI- 1.1*100. 

1LJ1  ( K,  M )  =  (  Q  (Kx^L/_XMJN(Hl-l  ,J*LQ£. _ 

WRITE! IOUT, 5  IK,! S! K, M I , PBLANKIK , M) ,H=1 ,101 
J  E-_!-X  L- 13  «_U2iJ3-x32_ 

32  SF!K)  =  SF!K)«-S.!K,1  I 
_ S A F l K )=SAF( KltClK.JLl 

33  CONTINUE 
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COMPUTER  PROGRAM  ORCBP  (FORTRAN  IV,  IBM  360/75) 

_ WRITE! I  OUT  §  9 ) _ 

WRITE (I  OUT, 1C)XM0N(I I ),N,NNAME 

_ DO  34  Kg  1 , KMAX _ 

WRITE ( IOUT,5)K,(0( K,M)fNBLANK(K,M) ,M=l,10) 

34  C ONTJ NUE _ 

DO  35  M= 1,10 

35  XMINIMI=9999999* _ _ 

C  CALCULATES  DEPARTURES  FROM  MEDIAN 

_ DO  45  K  =  3,KMAX,2 _ 

L=N-K+1 

_ _ 

J2=K. 

_ MID= I K+ 1 ) /2 _ 

DO  40  J J=  l ,  L 

_ J5= J2- 1  _ . _ 

DO  36  J= J 1 , J2 

36  Z  ( J  1  =  X(  LL  ,  J ) _ . _ 

DO  39  I  I  I = J 1 ,M1D 

_ DO  38  J=  J  1 ,  J  5 _ 

IF(ZIJ)-Zt J*l))38,38,37 
__  _3_7_H0LD=ZC  J ) _ 

zcji*zc  j+i'i” 

_  Z  ( J  +  l )=HOLD _ 

38  CONTINUE 

39  J  5= J  5- 1 _ 

Y(JJ)=Z(MID) 

_ J1=J1*1  _ _ 

J2=J2+1 

_ 40  MI  0=MID*JL _ 

DO  44  M=1 ,10 

_ SSUMgO> _ 

NMK 1=N— M-K+ 1 

_  DO  41  J=  1  ,NMK  1 _ _ 

JKMl=j+K*M-i 

41  SSUM=SSyMjfABSJ.YlJ)j1X(LL,JKMin 
C2=NMK1 

_ S I K , M )  =  SSUM/C2 _ 

IF (S(K,M)-XM  INI  Ml) 42,43,43 

_ 42  XMIN(M»  =  S(K,M^) _ 

L'CCIM  j  =  K 

_ 43  NBLANKI K, M) =NBL _ 

44  CONTINUE 

45  CONTINUE _ 

00  46  M=1 ,10 

_  K  K=  L  0  C  (  M  ) _ _ _ 

46  'NBLANKI  KK,M)  =NST 

WRITEUOUTV?)  _ _ 

. WRITE!  I  OUT,  14  )  XMONI I  IJ,N,NNA ME 

DO  47  K=3,KMAX,2 

47  WRITE! I0UT,5)K,I SI K,M) ,NBLANKIK,MI ,P=l,10l 
WR I TE 1 1  OUT , 9  ) 

WRTtEUOUT,  12)  XMONI  1 1  )  ,  N,NNAME 

DO  50  K  =  3,KMAXjl2 _ _ _ 

DO  48  M= 1,1 0 

48  S(K,M)~IS(K,MI/XMIN(M)-1. 1*100* 

.  WR II L 1 1  OUT , 5 ) K , l SI K , M ) , NBL ANK I K ,M I  ,  M=1 , 10 ) 

IF  1X1-13. )49, 50,49  _  _  _ 

49  s c f  i  k  i  =  sd"f  I* K  I  ♦  S I  K  ,  1 1 


FUg*  3 
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COMPUTER  rROCHUM  ORCBP _  (FORTRAN  IV,  IBM  360/75) 


5Q.XQN.TINUF. 

51  CONTINUE 

_ LI =13 _  __ 

WRITE ( I  OUT »  9 ) 

_ JbLRI  TE  ( I  OUT ,  2  1XK0NLII_J*N»_NNAM£. 

DO  52  K=1,KMAX,2 


53  CONTINUE 


STOP. 


r 


DATA  FUR  JAN 

OF  A  100 

YEAR  PERIOD 

SANTA  BARBARA,  CALIF. 

PRECIPITATION 

( INCHES) 

1868-1967 

3.97 

3.26 

0.25 

0.86 

2.53 

0.58  " 

”  4.54 

14.94 

7.56 

2.72 

7.17 

5.24 

1.30 

2.83 

1.13 

2.18 

6.3  3 

1  .23 

5.12 

oTTi 

10.15 

C .  29 

5.32 

0.45 

1.10 

4.41 

0.99 

6.25 

6.84 

4.35 

0.63 

4.48 

2.32 

4.d6 

1.36 

2.06 

0.46 

3.73 

4.26 

12.46 

4.29 

15.67 

2.91 

14.21 

C.4? 

3.14 

15.91 

4.94 

17.24 

3.25 

0.51 

1.20 

0.33 

5.  32 

4.64 

1.96 

1  .63 

9.60 

2.08 

1  .94 

0.0 

1.53 

5.82 

4.25 

2 . 4  C 

6.42 

1.49 

4.10 

0.73 

3.09 

1.90 

2.84 

6.39 

9.  b  8 

0.  BC 

12.84 

1  .44 

0.60 

0.40 

0’l60 

0.0 

1.40 

2.54 

2.83 

13.89 

1.78 

5.98 

4.  39 

7.19 

5.39 

3.71 

2.68 

3.12 

1.81 

2.  1  8 

1.79 

1.45 

0.76 

1.51 

7.61 

l 
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SIK, 

M1S0  FOR 

JAN  OF  A 

100  YfcAR 

PERIOD  SANTA  BARBARA 

,  CALIF. 

PRECIPITATION  ( INCHES  1 

1868-1967 

1 

2 

3 

4 

5 

6 

7 

8 

10 

K 

1 

29.95 

21.22 

30.80 

30.58 

34.35 

30.86 

32.09 

3".  19 

32.7  3 

2 

18.26 

18.50 

24.4? 

24.28 

24.78 

24.92 

23.06 

23.21 

2  3.50 

24.30 

3 

18.49 

19.39 

22.55 

23.62 

22.81 

22.36 

21 .23 

22.?  4 

22.08 

23. ?4 

4 

18.81 

18.86 

22.  14 

22.03 

20.49 

20.67 

20.62 

21  .00 

21.40 

21.84 

5 

18.33 

19.07 

2C.96 

19.73 

19.22 

20.11 

19.67 

20.45 

20.47 

21.34 

6 

18.60 

18.60 

18.97 

18.78 

19.01 

19.  39 

19.45 

19.86 

29.31 

?0. 6° 

7 

18.19 

16.85 

18.13 

18.60 

18.43 

19.  16 

19.03 

19.76 

19.77 

?C  ,55 

a 

16.61 

16.40 

18.10 

18.  19 

18.40 

18.87 

19.10 

o 

• 

19.8  > 

20.07 

9 

16.34 

16.54 

1  7.86 

18.25 

18.16 

18.94 

18.78 

19.52 

19.63 

?0 . 24 

10 

16.48 

16.52 

18.00 

18.05 

18.32 

18.67 

18.97 

19.37 

19.82 

19.68 

u 

16.46 

16.74 

17.85 

18.21 

18.12 

18.86 

18.87 

19.58 

19.34 

19.73 

12 

16.68 

16.68 

18.01 

16.02 

18.34 

18.81 

19.18 

19.39 

18.77 

13 

16.61 

16.90 

17.84 

19.23 

18.33 

19.01 

18.76 

19.21 

18.51 

18.88 

14 

16.85 

16.84 

18.08 

18.24 

18.57 

18.72 

18.84 

18.36 

18.65 

18.03 

15 

16.79 

17. 12 

16.  T9 

16.47 

18.33 

18.80 

18.03 

18.50 

18.76 

18.90 

16 

17.  C8 

17.18 

16.34 

18.26 

18.45 

18.00 

18.20 

18.61 

M.71 

18.96 

17 

17.13 

17.46 

18.  14 

18.37 

17.66 

18.15 

18.31 

18.58 

ld.7« 

18.82 

16 

17.39 

17.33 

18.  ?4 

17.64 

17.82 

18.28 

18.28 

18.65 

18.67 

18.71 

19 

17.27 

17.47 

17.81 

17.79 

17.94 

18.26 

18.37 

18.54 

13.58 

18.74 

20 

17.41 

16.81 

17.67 

17.94 

17.95 

18.36 

18.30 

18.48 

19.64 

18. 8H 

21 

16.74 

17.00 

17.81 

17.96 

18.06 

18.30 

18.27 

18.54 

18.76 

18.93 

22 

16.92 

17.19 

17.82 

18.07 

18.05 

18.27 

18.37 

18.70 

19.85 

19.0b 

23 

17.10 

17.23 

17.95 

18.06 

18.04 

18.36 

18.52 

18.79 

18.99 

1 0  •  c  \ 

24 

17.14 

17.37 

17.95 

18.05 

18.15 

18.51 

18.63 

18.92 

13.9b 

19.20 

25 

17.29 

17.40 

17.96 

18.15 

18.29 

18.62 

18.77 

18.91 

19.15 

19. 2" 

26 

17.32 

17.43 

18.06 

18.28 

18.41 

18.78 

18.77 

19.11 

19.17 

19.4" 

27 

17.37 

17.55 

18.20 

18.41 

18.57 

18.  79 

18.98 

19.13 

19.37 

19.33 

28 

17.48 

17.70 

18.34 

18.58 

18.59 

19.00 

19.01 

1«.34 

19.31 

19.89 

29 

17.64 

17.85 

18.50 

16.61 

18.  R1 

19.05 

19.22 

19.30 

19.56 

19.71 

30 

17.79 

18.01 

18.54 

18.83 

18.87 

19.25 

19.18 

19.53 

19.69 

20.05 

31 

17.95 

18.07 

16.  76 

18.90 

19.07 

19.2! 

19.41 

19.67 

20.0  3 

19.05 

32 

18.01 

18.29 

18.83 

19. n 

19.06 

19.45 

19.55 

20.01 

19.  op 

19.15 

33 

18.25 

18.39 

19.04 

19.10 

19.29 

19. 60 

19.90 

18.99 

19.11 

16.97 

34 

18.34 

18.60 

19.04 

19.33 

19.46 

19.95 

18.97 

19.11 

16.92 

17.17 

35 

18.56 

18.63 

19.28 

19. 5J 

19.80 

18.94 

19.  0" 

16.93 

17.13 

15.67 

36 

18.59 

16.87 

19.45 

19.85 

18.78 

19.07 

16.80 

17.14 

15.52 

15.61 

37 

18.84 

19.06 

19.60 

18.82 

18.92 

1  6  •  36 

17.02 

15.50 

15.45 

15.78 

38 

19.03 

19.40 

18.74 

16.96 

16.67 

17.07 

15.35 

19.47 

15.69 

13.0? 

39 

19.38 

18.36 

18.89 

16.71 

16.89 

15.40 

15.3? 

15.71 

12.98 

13. ?P 

HI 

18.34 

18.51 

16.66 

It. 94 

15.23 

15.37 

15.5b 

13.01 

13.23 

9.96* 

41 

18.50 

16.37 

16.87 

15.35 

15.18 

15.61 

12.96 

13.27 

10.02* 

lr  •  ?9 

42 

16.37 

16.59 

15.30 

15.30 

15.42 

13.02 

13.23 

10.07* 

1".  36 

10.29 

43 

16.58 

15.12 

15.22 

15.54 

12.96 

13.30 

10.17* 

10.39 

10.32 

10.46 

44 

15.10 

15.05 

15.46 

13.07 

13.24 

10.22* 

10.49 

1C. 35 

1  2.49 

10.33 

45 

15.00 

15.28 

13.09 

13.35 

10.28* 

10.54 

10.42 

10.45 

l*.  32 

10.65 

46 

15.23 

12.88 

13.37 

1C. 35* 

10.59 

10.48 

10.57 

10.34 

10.63 

10.76 

47 

12.83 

13.16 

10. 35* 

1C. 66 

10.53 

10.63 

10.42 

13.75 

10.87 

48 

13.11 

10.26* 

10.65 

1C. 57 

10.65 

10.44 

10.74 

10.78 

10.83 

10.89 

49 

10.22* 

10.56 

10.56 

10.68 

10.47 

10.76 

r*- 

O 

10.85 

10.81 

_5fl_ 

10.50 

_ 10.42 

10.63 

1C. 45 

10.79 

10.90 

10.93 

Mtl ,  • 

10.73 

10.87 

90 
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r 


F(K,K) 

FOB  JAN  OF  A  ICO  YEAR 

PERIUO 

SANTA  BARBARA, 

CALIF. 

PRECIPITATION  (INCHCS) 

1868-1967 

1  2  3 

4 

5 

6 

7  8  9 

10 

K 


1 

192. 9A 

IC6.72 

197.59 

215.55 

197.e3 

235. 98 

203.59 

218.84 

2011 .25 

228.78 

2 

78.66 

60.29 

135. 5A 

13A.71 

141.21 

143.79 

126. BO 

1  30 . 6  3 

134.51 

144.09 

3 

80. 8A 

68.90 

117.87 

128.2b 

121.98 

118. V  5 

lCb.dO 

121.01 

120.31 

133.37 

A 

8A.02 

83.77 

113.HR 

112.96 

99.37 

102.15 

102.79 

1  OR . 62 

113.58 

119.31 

5 

79.28 

85  •  85 

102.53 

90.73 

b7.06 

96.66 

93.51 

K3.22 

l°4.2l 

1  14. 36 

6 

81.99 

81.22 

83..:  7 

81. 50 

85,00 

89.65 

91.32 

97.27 

102.66 

106.90 

7 

77.96 

6A.17 

75.19 

79.82 

79.3b 

87.42 

87.19 

96.29 

,’.26 

106.45 

8 

62. 5C 

59.76 

7a.  69 

75.79 

79.03 

84  •  ft7 

8  7.83 

92.64 

93.44 

If  1 .63 

9 

59. 8A 

61.13 

72.79 

76.35 

76.78 

85.29 

84.70 

93.90 

9ft. 84 

103.32 

10 

61. 2A 

61.01 

73.93 

74.47 

78.26 

82.63 

Bb  .  5o 

92.43 

97.81 

07.62 

11 

61.00 

63.12 

72. Aft 

75.97 

76.39 

64 . 52 

35.66 

94  .  SO 

93.03 

98.14 

12 

63.13 

62.52 

7A.01 

7A.20 

7b. 51 

84.00 

87.78 

90.52 

93.47 

88. 50 

13 

62.52 

64 .  7 1 

72. AC 

7C.22 

78. AA 

65.91 

8A.52 

90.87 

84.68 

89.60 

l  A 

6A.8C 

6A.09 

7  A ,  6  3 

76.27 

d0.76 

83. 10 

85.34 

82.40 

86.10 

90.59 

15 

6A.29 

66.83 

7A.  77 

7 1  •  5A 

78.42 

8  3.91 

77.35 

83.81 

87.16 

89.83 

_Ul 

67.06 

67. A2 

77.  15 

76.  A6 

79.54 

76.  93 

78.99 

84.91 

36.74 

90.45 

17 

67.57 

70. 1A 

75.26 

77.58 

71.90 

77.57 

90.09 

84.58 

37.42 

88.  OR 

18 

70.17 

68.88 

76.26 

70.50 

73.40 

78.  7B 

79.83 

85.25 

86.27 

87.96 

19 

68.96 

70.21 

69.  1  6 

71.9ft 

74.61 

78.59 

80.65 

34.22 

85.44 

tP  .  25 

20 

70.3b 

63.  79 

70.70 

7  3. 3o 

74.67 

79.53 

80.04 

8  3.65 

86.04 

89.53 

21 

63.79 

65. 6A 

72. OA 

73.55 

75.50 

78.96 

79.73 

84.17 

87.20 

90.15 

?? 

65.50 

67. AS 

72.  19 

7  4  •  6  3 

75.69 

78.71 

00.73 

36.76 

88.12 

91.39 

23 

67.26 

67.92 

73.  38 

74.54 

75.61 

79.ftb 

82.19 

86.71 

89.47 

90.95 

_2A 

67.68 

69.25 

73.  AA 

74.45 

76.43 

81.  JO 

83.23 

88.  p7 

89.10 

92.81 

25 

69.11 

69.51 

73.55 

75.41 

78. "A 

M2.  15 

84  .67 

87.63 

91.12 

92. 86 

26 

69.  AA 

69.87 

7a.  Ab 

76.72 

79.19 

S3. 66 

84 . 4‘3 

89.82 

91.24 

«4.8  2 

27 

69.88 

71.01 

75.  65 

7  7.  97 

80.7  6 

63.  75 

d6  *6  7 

90.11 

93.23 

94.10 

26 

71.00 

72.  AA 

77.15 

75.61 

80.67 

35.  86 

87.  CA 

92.1b 

92.71 

96. 73 

29 

72.52 

73.92 

76.73 

79.92 

83.04 

84.  13 

89. C6 

91.73 

95.1" 

98  .fo 

74,05 

75.  A8 

70.13 

81.96 

83.o7 

88. 24 

88.65 

94.04 

96.49 

10.39 

31 

75.61 

76.06 

ei.2i 

82.63 

d5.61 

ft7.  -14 

60.42 

05.42 

99.89 

91.31 

32 

76.23 

78.27 

61. 9C 

84.61 

85.4'j 

90.24 

42.35 

98.76 

39.60 

92.31 

33 

78.50 

79.15 

63.91 

84.57 

87.77 

91.76 

05.73 

8b.  67 

90.72 

7C.45 

3A 

79. A3 

81.27 

83.95 

86.88 

89.35 

95.10 

•  6  *) 

89.84 

63.85 

72.44 

35 

81.57 

81. A9 

86.27 

8e.52 

62.7  3 

85.2  3 

66.93 

68.25 

70.89 

56. 3P 

36 

81.8b 

83.86 

87.95 

9  1 . 6*. 

82.75 

C6.52 

66.30 

70.31 

5<»  •  66 

55.78 

37 

8A.35 

85.70 

91.31 

81.87 

84.14 

64  •  'IV 

67.43 

54.03 

54.22 

58.20 

38 

86.17 

H9.07 

61.09 

b  3 . 2  3 

6?  •  2 b 

6  7.  )0 

50.66 

53.6? 

56.61 

3f  .  84 

39 

89.58 

78.85 

d2 . 52 

6  1.56 

64.42 

50.64 

50.6  7 

56.06 

29.52 

33.42 

AO 

79.37 

80.36 

6C.53 

6  3.69 

48.19 

50.  34 

53.^6 

29,?1 

32.06 

O.C  * 

A 1 

80.98 

59.50 

63. OA 

4  8  .  ft  o 

47.7'' 

52.70 

27.48 

31.3b 

o.O  * 

3.39 

A2 

60.12 

61.69 

A  7. 7  7 

47.91 

5  C  .09 

27.  IB 

'  .  1  3 

* 

3.36 

3. 36 

A3 

62. 2A 

A7.29 

A7.C1 

5  C  .  2  0 

26.15 

30.  r 

0  * 

3.24 

2.93 

5.08 

AA 

AT. 68 

A6.6C 

A9. 34 

2  t  .  3  A 

28.85 

1  ♦ 

3.1 ; 

2 .  Pt) 

4.64 

3.72 

A5 

A6.70 

A8.86 

26.  A7 

2  5 . 06 

C.O  ♦ 

?.  34 

2.46 

4 . 24 

2.95 

6.93 

A6 

A9.03 

25.  A8 

29.  1H 

c.o  * 

3.05 

2.47 

?  .94 

2.73 

6.12 

8.10 

A7 

25.  A8 

28.21 

0.0  • 

3.01 

2.51 

3.93 

2.49 

5.14 

7.29 

9.16 

A8 

28.29 

0.0  * 

2.90 

2.12 

3.66 

2.13 

6  .ft" 

7.11 

3.12 

9.36 

A9 

O.C  * 

2.89 

1.90 

3.22 

1.85 

6.38 

0 . 9  ft 

7.92 

8.31 

8.57 

50 

2 • 7A 

1.5A 

2.72 

C.09 

5.04 

O  • 

7.4  7 

7.71 

7.03 

9.21 
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Q  IK.IO 

FOR  JAN  OF  A  ICO  YEAR  PERIOD  SANTA  r.AKBARAf  CALIF. 

PRECIPITATION  1  INCHES!  1866- 

M>  1 

2  3  4  5  6 

7  8  9  1  r 

1 

3.93 

3.02 

4. 02 

4.10 

3.98 

4.  35 

4.01 

4.25 

3.8  3 

4  •  1  3 

2 

3.09 

3.  19 

3.7b 

3.74 

3.83 

3.94 

3.77 

3.72 

4.  7° 

3.74 

3 

3.18 

3.24 

3.5b 

3.70 

3.72 

3.  81 

3.52 

3.63 

3.40 

3.64 

A 

3.23 

3.18 

3.55 

3.62 

3.60 

3.59 

3.54 

3.47 

4  •  4h 

3.52 

5 

3.20 

3.24 

3.50 

3.53 

3.45 

3.55 

3.39 

3.44 

j  .  ®b 

3.46 

6 

3. 23 

3.25 

3.41 

3. 38 

3.45 

3.40 

3.39 

3.3  7 

3.  36 

3.3/ 

7 

3.21 

3.2C 

3.30 

3.40 

3.33 

3.  39 

3.35 

3.35 

3.28 

3.  3b 

8 

3.15 

3.  13 

3.  33 

3.24 

3.34 

3.34 

3.35 

3.26 

3.2« 

3.23 

9 

3.11 

3.  16 

3.25 

3.32 

3.29 

3.33 

3.24 

3.26 

3.21 

3.22 

10 

3.15 

3.  12 

3.29 

3.28 

3.29 

3.24 

3.26 

3.2  4 

3.21 

3.  14 

11 

3.11 

3.  15 

3.27. 

3.29 

3.21 

3.25 

3.22 

3.22 

3.14 

3.  14 

12 

3.15 

3.14 

3.28 

3.21 

3.22 

3.22 

3.22 

3.1  7 

3.13 

3.06 

13 

3.13 

3.  16 

3.20 

3.21 

3.20 

3.22 

3.18 

3.  16 

4.04 

3.'  4 

1* 

3.15 

3.10 

3.21 

3.18 

3.21 

3.18 

3.17 

3.07 

>. 05 

3.07 

15 

3.11 

3.12 

3.18 

3.19 

3.17 

3.17 

3.07 

3."7 

3.08 

3.04 

16 

3.13 

3.  10 

3.  19 

3.16 

3.17 

3.07 

3.07 

3.11 

3.0b 

3.04 

17 

3.11 

3.12 

3.17 

3.16 

3.06 

3.06 

3.10 

3.07 

OS 

3.04 

18 

3*12 

3.  10 

3.  16 

3.06 

3.05 

3.  n 

3.06 

3.0/ 

3.06 

3.01 

19 

3.10 

3.09 

3. 05 

3.05 

3.08 

3.06 

3.06 

3  .r  7 

3 . 0  3 

3.05 

20 

3.10 

3.00 

3.C5 

3.08 

3.05 

3.06 

3.07 

3.09 

3.0b 

3.06 

21 

3.00 

3.00 

3.07 

3.05 

3.C6 

3.0b 

3.05 

3.0  9 

3.07 

3.19 

22 

2.99 

3.03 

3. 04 

3.05 

3.08 

3.  06 

3.09 

3.°9 

3.12 

.3.  "6 

23 

3.02 

3.00 

3.05 

3.07 

3.0b 

3.09 

3.10 

3.14 

3.  i  i 

3.ii 

2* 

2.98 

3.01 

3.07 

3.C6 

3.09 

3.  1 

3.14 

3.13 

3.13 

3.12 

25 

2.99 

3.03 

3.05 

3.08 

3.09 

3.14 

3.13 

3.14 

3.15 

3.13 

2b 

3.01 

3.02 

3.08 

3.98 

3.12 

3.13 

3.14 

3.15 

3.15 
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